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1 Introduction 

The primordial inflation [1—3] has two key predictions - creating the scalar density per¬ 
turbations and the tensor perturbations during the accelerated phase of expansion [4], 
for a review, see [5]. One of the predictions, namely the temperature anisotropy due 
to the scalar density fluctuations has now been tested very accurately by the obser¬ 
vations from the temperature anisotropy and polarization in the cosmic microwave 
background (CMB) radiation [6-8]. In the last year the detection of tensor modes 
has been initially confirmed by the ground based BICEP2 experiment [9], which ac¬ 
cording to the initial claim has detected for the first time a non-zero value of the 
tensor-to-scalar ratio at 7cr C.L. The value obtained by the BICEP2 team in con¬ 
junction with Planck (2014) 1 WMAP-9 I high L+BICEP2 (dust) 1 to put a bound 

throughout the article we use red for the obsoleted Planck (2014)+WMAP-9+high L+BICEP2 
(dust) data. 
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on the primordial gravitational waves, via tensor-to-scalar ratio, within a window: 

0.15 < r{K) = P T {h)/P s {h) < 0.27, (1.1) 

at the pivot scale, fc* = 0.002Mpc _1 [9], where Pt and Ps denote the power spectrum 
for the tensor and scalar modes, respectively. But just after releasing this result BI- 
CEP2 analysis was put into question by several works [10-13] on its correctness in 
the physical ground. Most importantly it accounting for the contribution of fore¬ 
ground dust which will shift the value of tensor-to-scalar ratio r downward by an 
amount and further better constrained by the joint analysis performed by Planck 
and B1CEP2/Keck Array team [14]. The final result is expressed as a likelihood 
curve for r, and yields an upper limit: 

r(K) = P T (K)/P s (h)< 0.12 (1.2) 

at the pivot scale, k * = 0.05Mpc _1 with 2 a confidence. Marginalizing over dust 
contribution and r, finally it is reported that the lensing B-modes are detected at 
la significance. Very recently in ref. [15] the Planck team in 2015 data release also 
fixed the upper bound on the tensor-to-scalar ratio as: 

r(h) = P T (h)/P s (h)< 0.11 (1.3) 

at the pivot scale, fc* = 0.002Mpc -1 with 2cr C.L. and perfectly consistent with the 
joint analysis performed by Planck and BICEP2/Keck Array team. 

Note that large r(/c*) is possible if the initial conditions for gravitational waves is 
quantum Bunch-Davis vacuum [16] 2 , for a classical initial condition the amplitude 
of the gravitational waves would be very tiny and undetectable, therefore this can 
be treated as the first observable proof of quantum gravity. However, apart from the 
importance and applicability of quantum Bunch-Davis vacuum on its theoretical and 
observational ground it is still not at all clear from the previous works in this area 
that whether the quantum Bunch-Davies vacuum is the only source of generating 
large value of r(fc*) during inflationary epoch or not. One of the prime possibilities 
comes from the deviation from quantum Bunch-Davies vacuum aka consideration 
of quantum non-Bunch-Davies or arbitrary vacuum in this picture which may also 
responsible for the generation of large r(fc*) during inflation 3 . 

In this paper, our aim will be to illustrate that it is possible to explain the current 
data sets within a sub-Planckian model of inflation, where: 

• (j) o M p - vev of the inflaton must be bounded by the cut-off of the particle 
theory, where M p = 2.4 x 10 18 GeV. We are assuming that 4 dimensions M p 
puts a natural cut-off here for any physics beyond the Standard Model. 

2 Apart from the correctness of this argument, it is additionally importnat to mention here that, 
still it doesn’t require the initial conditions for the quantum vacuum to be Bunch-Davis strictly. 

3 In this article I have not explored the possibility of non Bunch Davies vacuum. 
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\A(f)\ Ri |0* — 0 e | < M p - the inflaton potential has to be flat enough during 
which a successful inflation can occur. Here 


0* > 00 > 0e (1-4) 

represents the field VEV, and A0 denotes the range of the field values around 
which all the relevant inflation occurs, 0* corresponds to the pivot scale and 0 e 
denotes the end of inflation. Note that the flatness of the potential has to be 
fine tuned and also there is no particle physics symmetry which can maintain 
the flatness. We will assume 

W'(0o) « 0 , (1.5) 

where V(0) denotes the inflaton potential, and prime denotes derivative w.r.t. 
the 0 field. Naturally, the potential has to be flat enough within A0 to support 
slow roll inflation. 


The above requirements are important if the origin of the inflaton has to be 
embedded within a particle theory, where inflaton is part of a visible sector gauge 
group, i.e. Standard Model gauge group, instead of an arbitrary gauge singlet. If the 
inflaton is gauged under some gauge group, as in the case of a minimal supersym¬ 
metric Standard Model (MSSM), Ref. [17], then the inflaton VEV must be bounded 
by Mp, in order to keep the sanctity of an effective field theory description 4 . 

Our prescribed methodology will be very generic, with a Taylor expanded po¬ 
tential around VEV 0o is given by 5 : 


Tr( ,\ (0 ~ 4>o) n fd n V(4>)\ 


n =0 


= V(M + V'(0o)(0 - 00 ) + V^(0 - 0o) 2 + AAV - 0o) 3 

C""(0o) , 

H-——(0 - 0o) H- j (1-6) 

where the expansion co-efficients are characterized as follows: 


• The first term: 

VM « Mp 4 (1.7) 

denotes the height of the potential. Also this term will play the most significant 
role in fixing the scale of inflation within the present framework. 

4 An arbitrary moduli or a gauge singlet inflaton can take large VEVs ( super-Planckian ) as 
in the case of a chaotic inflation [2]. Although, in the case of assisted inflation [18], see chaotic 
assisted inflation [19], the individual VEVs of the inflatons are sub-Planckian. 

5 For the sake of completeness I suggest the readers to see ref. [20] , where I have explicitly studied 
the inflationary reconstruction technique within the framework of Randall Sundrum single brane 
set up, using the observational constraints obtained from Planck 2015 and BICPE2/Keck Array 
joint constraints. 
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The Taylor expansion coefficients of the effective potential: 



<K 3 , 

(1.8) 

V"(ct> o) 

<M 2 p , 

(1.9) 

V m (fa) 

< M p , 

(1.10) 


< 0(1), 

(1.11) 


determine the shape of the potential in terms of the model parameters. 

• The prime denotes the derivative w.r.t. <fi. In particular, some specific choices 
of the potential would be a saddle point , when 

H0 O ) = 0 = l/"(0o), (1-12) 

an inflection point , when 

V"^) = 0. (1.13) 

Previous studies regarding obtaining large r(kflj within sub-Planckian VEV mod¬ 
els of inflation have been studied in Refs. [21-23]. In Ref. [21], the authors could 
match the amplitude of the power spectrum, Pg, at the pivot point, but not at the 
entire range of A0 for the observable window of AN, where N is the number of 
e-foldings of inflation. In Ref. [24], the authors have looked into higher order slow 
roll corrections by expanding the potential around 0 O . They pointed out that large 

r ~ 0.05 (1.14) 

could be obtained in an inflection-point model of inflation where the slow roll param¬ 
eter, ey, changes non-monotonically (for a definition of ey, see Eq. (2.8)). The ey 
parameter first increases within the observational window of AN and then decreases 
before increasing to exit the slow roll inflation by violating the slow roll condition, 
i.e. ey ~ 1 6 . The value of r was still small in order to accommodate the WMAP 

6 The smallness of the observational window on AN also allows us to exploit another loophole 
in the derivation of the Lyth bound [25], namely the assumption that ey increases monotonically. 
This seems a natural assumption given that during inflation 

ey < 1 (1.15) 

and at the end of inflation 

ey ~ 1, (1.16) 

current observational constraints also strongly favour ey increasing during the AN ~ 17 e-folds 
of the window achieved by the CMB distortion observation. However, outside this window the 
behaviour of ey is not constrained and the assumption of monotonicity is not strictly necessary. 
Relaxing this assumption makes it possible to construct a scalar potential for a single field that 
violates the Lyth bound. For more details see Ref. [24], 
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data, which had probed roughly AN fa 8 as compared to the Planck, which has now 
probed AN fa 17 e-foldings of inflation [23, 26-29] can be achieved by the CMB dis¬ 
tortion observation. A generic bound on tensor-to-scalar ratio was explicitly derived 
in Ref. [23, 26, 27] for a generic type of inflection-point inflationary model, where I 
set R"(0o) = 0 in the Taylor expanded form of the potential. 

In this paper I will consider the full potential of Eq. (1.6), and onr prime objec¬ 
tive is to determine the values of the Taylor expansion co-efficients V(4>o), V'(4>o), 
V"((j) 0 ), V'"((f) o) and 0 ) from the latest Planck 2015 and BICEP2/Keck Ar- 
ray+Planck 2015 joint data 7 . In this respect I will be reconstructing the inflationary 
potential around 0 o and at the pivot scale, 0* = 0(fc*), where I fix £;* = 0.002Mpc _1 . 
We will also provide for the second order consistency relations for a sub-Planckian 
excursion of the inflaton field. This will be treated as an observational discriminator 
which could rule out various sub-Planckian models of inflation in future. 

Within the region of AN fa Oifl — 17) e-foldings [23, 26-29], I will be able to 
constrain the power spectrum: P$, spectral tilt: ns, running of the spectral tilt: as, 
and running of running of the spectral tilt: Ks, in the background of ACDM model 
for: 

Planck (2013)+WMAP-9+high L data sets: [7, 8] 


r(k it ) < 0.12 (within 2cr C.L .), (1.17) 

ln(10 10 P 5 ) = 3.089lolo27 (within 2d C.L.), (1.18) 

ns = 0.9600 ± 0.0071 (within 3a C.L), (1.19) 

as = clns/dlnk = —0.013 ± 0.009 (within 1.5a C.L), (1-20) 

Ks = d 2 ns/d In k 2 = 0.020lo!oi5 (within 1.5a C.L.). (1.21) 


Planck (2014)+WMAP-9+high L+BICEP2 (dust) data sets:[9] 


0.15 < r(kfl) < 0.27 (1.22) 

ln(10 10 P 5 ) = 3.O890 o|o 27 (within 2a C.L), (1.23) 

n s = 0.9600 ± 0.0071 (within 3a C.L), (1.24) 

as = dns/dlnk = —0.022 ± 0.010 (within 1.5a C.L), (1-25) 

Ks = d 2 ns/d In k 2 = 0.020lo!oi5 (within 1.5a C.L.). (1.26) 


' Additionally I have also mentioned the results using the obsoleted Planck+WAMP-9+high 
L+BICEP2(dust) data throughout the paper to explicitly show the validity of our prescribed 
methodology for any kind of observed data set. 
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Planck (2015)+WMAP-9+high L(TT) data sets:[15] 

r(k i ) < 0.11 (within 2 a C.L. ), (1-27) 

ln(10 10 P 5 ) = 3.089 ± 0.036 (within 2a C.L), (1.28) 

n s = 0.9569 ± 0.0077 (within 3a C.L), (1.29) 

as = dns/d In k = O.OII^oq^ (within 1.5a C.L), (1.30) 

Ks = d 2 ns/d In k 2 = 0.029^o;o]q (within 1.5a C.L) . (1-31) 

Planck (2015)+BICEP2/Keck Array joint data sets:[14] 

r{k±) < 0.12 (within 2a C.L), (1.32) 

ln(10 10 P 5 ) = 3.089^027 (within 2a C.L), (1.33) 

n s = 0.9600 ± 0.0071 (within 3a C.L), (1.34) 

as = dns/dlnk = —0.022 ± 0.010 (within 1.5a C.L), (1.35) 

Ks = d 2 ns/d\nk 2 = O.C^O^o}® (within 1.5a C.L) . (1.36) 


In this paper- 

• We will briefly recap the key equations for the inflationary tensor-to-scalar ratio 
in the most generalized case by taking into account of the effect of higher order 
slow-roll corrections. Throughout the analysis of the paper I assume: 

1. Inflaton field </> is minimally coupled to the Einstein gravity sector. 

2. Slow-roll prescription perfectly holds good after considering higher order 
slow-roll corrections. 

3. Convergence of the Taylor expanded potential. 

4. For the numerical estimations the generic form of Taylor expanded poten¬ 
tial is truncated at the fourth power of inflaton field. 

5. We also assume that the major contribution for the generic version of 
inflationary potential comes from the first term V((j) o) of the Taylor series 
for which convergence of the taylor sereies holds good perfectly within the 
present framework. We will explicitly show in the next sections, within 
the slow-roll regime of inflation this assumption is compatible with the 
results obtained by applying the recent observational constraints. 

6. As the contribution from all the non-renormalizable effective field theory 
operators are highly suppressed by the various powers of the UV cut-off 
scale of the effective field theory Ajjv, in the present context I neglect all 
such contribution due to its smallness. 

7. Initial condition for inflation is fixed via Bunch-Davies vacuum. 
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8. Sound speed is fixed at cs = 1. 

9. UV cut-off of the effective theory is fixed at A uv = M p , where M p is the 
reduced Planck mass. But in principle one can fix the scale between GUT 
scale and reduced Planck scale i.e. Aqut < A uv < M p . But in such a 
situation A uv acts as a regulating parameter in the effective field theory. 
To avoid all such complications I fix it at M p . 

• We will then derive the most general bound on r(k ir ) for a generic sub-Planckian 
VEV inflation by considering the effect of running and running of the running 
in primordial scalar and tensor power spectrum, and the corresponding values 
of H+ and V (0*) within the framework of effective field theory. For complete¬ 
ness of the presented analysis, in the Appendix of this paper I give also the 
expression for the bound obtained from the various parameterization in the 
primordial power spectrum. 

• We will then discuss the non-monotonic behaviour of slow-roll parameter and 
its effectiveness to evade the Lyth bound [25] in section 4. 

• Further I reconstruct the shape of the potential in section 5, by providing the 
observational constraints from Planck (20 13 ) : W.YI AIMI—high- 1 . [7, 8], Planck 
(2014)+WMAP-9+high-L+BICEP2 (dust) [9], Planck (2015)-f-WMAP-9 | high- 
L(TT) [15] and Planck (2015)+BICEP2/Keck Array joint constraints [14] on 
the various Taylor expansion co-efficients V'(0 0 )> V"(0 O ), V w (0o) and V' w (0o) 
at VEV 0o ; which are expressed in terms of the Taylor expansion co-efficients 
V'(0*), V"(0*), V"'(0*) and V""(0*) at pivot/CMB scale 0*. 

• In section 6, I will discuss the inflationary consistency relationships considering 
up to second order correction in slow-roll parameters. 

• In section 7, I will consider a specific case of inflection point inflation for the 
purpose of illustration. 

• In section 8, I will consider a specific case of saddle point inflation for the 
purpose of illustration. 

• Finally in section 9 I have discussed the multipole scanning of CMB TT, TE, 
EE and BB spectra via reconstructed potential by applying the constraints from 
Planck (2013)+WMAP-9+high-L, Planck (2014)+WMAP-9+high-L+BICEP2 
(dust), Planck (2015)+WMAP-9+high-L(TT) and Planck (2015)+BICEP2/Keck 
Array data. 

• Additionally, in the Appendix of our paper for completeness I will provide 
all the key equations and detailed derivations of main results which are used 
frequently throughout the prescribed analysis of the paper. 
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2 Brief introduction to tensor to scalar ratio in inflation 


The tensor to scalar ratio can be defined by taking into account of the higher order 
corrections, see Refs. [23, 30, 31]: 


r = lgej; - 11 - (Cg +1)^] 2 —_ _ 

[1 — (2 C E + 1 )ch + CeVh] 


( 2 . 1 ) 


where 

C E = 4(1h 2 + 7e )-5 (2.2) 

with 7 e = 0.5772 is the Euler-Mascheroni constant [30]. In Eq (2.1) the Hubble slow 
roll parameters (e#,?-/#) are defined as: 


din H 

H 

(2.3) 

d In a 

H 2 ’ 

d In 0 

0 

(2.4) 

Vh - „ 

am a 

770 ’ 


where dot denotes time derivative with respect to the physical time. Now considering 
the effect from the leading order dominant contributions from the slow-roll param¬ 
eters, the Hubble slow-roll parameters can be expressed in terms of the potential 
dependent slow-roll parameters, (ey, rjv), as: 


e H & e v -\ -, (2.5) 

Vh ~ Vv ~ e v H-, (2.6) 


where • • • comes from the higher order contributions of (ey, 7/y). 

The tensor to scalar ratio can be re-expressed in terms of inflationary potential 
as: 

r * 16e v - 11 - (Ce + hffi! 2 (2.7) 

[1 — (3 Ce + l)ey + Cet/v] 

where slow-roll parameters (ey, ?/y) are given by in terms of the inflationary potential 
V((p), which can be expressed as: 



r] V = Mp 



( 2 . 8 ) 

(2.9) 


We would also require two other slow-roll parameters, (£y, cry), in our analysis, which 
are given by: 


£ 


2 

V 


a 


3 

V 


M 4 

M 6 

UVJ-p 



V' 2 V""\ 
V 3 ) ' 


( 2 . 10 ) 

( 2 . 11 ) 


-8 - 









Note that I have neglected the contributions from the higher order slow-roll terms, 
as they are sub-dominant at the leading order. With the help of 


d _ _ ^ d 
dink p 1 — en d(j) 


- M r 


y/2€y d 


P 1 — ey d(j) ’ 

1 can derive a simple expression for the tensor-to-scalar ratio, r, as: 


r = 


8 (1 — e v ) 2 [1 — (Ce + 1 )tvY 


( d<t> V 

\dhikj ' 


M 2 [1 — (3 Ce + l)ey + CeVvY 
We can now derive a bound on r(k) in terms of the momentum scale: 


( 2 . 12 ) 


(2.13) 



dk r[k) 

TV TT~ 


= J_ f 4 ’* 11 (! ~ e v) t 1 ~ (£e + l)e v \ 

M P J<f, e [1 — (3 Ce + l)ey + C E r)v\ ’ 

— [ #(1 — e v ) [1 +C E (2e v - r)v) + ■■■■], 


M, 


P -He 


M p \ A (f) 


(2 C E — 1) / d(j) ey — C E / dcf) rjy 



(2.14) 


where note that 

A0^0*-0 e >O (2.15) 

is positive in Eq. (2.14), and (p e denotes the inflaton VEV at the end of inflation, 
and 0* denote the held VEV when the corresponding mode fc* is leaving the Hubble 
patch during inflation. Here 1 have used the slow-roll approximation 


<p/H ~ 'Jpy. 


(2.16) 


The physical significance of the Eq (2.14) are appended below: 

• This gives the the analytical expression for the held excursion during inflation 
in terms of tensor-to-scalar ratio and other inhationary observables. 

• This relation can be treated as a discriminator between sub-Planckian and 
super-Planckian inhationary models, depending on the value of held excursion. 

• This relation also justihes the validity and correctness of Effective Field Theory 
framework, depending on the value of held excursion. 

• Also this relation can be used to break the degeneracy between various class 
of inhationary models. 

8 We have derived some of the key expressions in an Appendix, see for instance, Eq. (10.8), 
which I would require to derive the above expression, Eq. (2.13). 
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Note that Acp > 0 implies that the left hand side of the integration over momen¬ 
tum within an interval, k e < k < k+, is also positive, where k e represents correspond¬ 
ing momentum scale at the end of inflation. Individual integrals involving ey and r/y 
are estimated in an Appendix, see Eqs. (10.12) and (10.13). 

In order to perform the momentum integration in the left hand side of Eq (2.14), 
I have used the running of r(k), which can be expressed as: 


where 


r(k) = r(k+) 


^ k ^ a +l ln (^-)+f ln2 (^-)+-- 


(2.17) 


a = riT — ns + 1, (2-18) 

b = (a T — as ), (2-19) 

c = (kt — ks) (2.20) 


defined at the momentum pivot scale fc*. These parameterization characterizes the 
spectral indices, ns, nr , running of the spectral indices, as, ax, and running of the 
running of the spectral indices, Ks, Kt- Here the subscripts, (S', T), represent the 
scalar and tensor modes. Now substituting the explicit form of the potential stated 
in Eq. (7.27), I can evaluate the crucial integrals of the first and second slow-roll 
parameters (ey, r/y) appearing in the right hand side of Eq. (2.14). For the details 
of the computation, see appendix. 

It was earlier confirmed by the WMAP9 | high-/1 BAO | 77 0 combined constraints 
that [6]: 

a s = -0.023 ±0.011, (2.21) 

k s = 0 ( 2 . 22 ) 

within less than la C.L. . After the Planck release it is important to see the impact 
on r(/c*) due to running, and running of the running of the spectral tilt by modifying 
the generic power law form of the parameterization of tensor-to-scalar ratio. The 
combined Planck (2013)+WMAP-9 constraint confirms that [8]: 

= -0.0134 ± 0.0090, (2.23) 

ks = 0.020±°;™j (2.24) 

within 1.5a statistical accuracy, which additionally includes 


k s ^ 0 (2.25) 

possibility for the first time. Also the recent combined Planck (2015)+WMAP-9 
constraint confirms that [15]: 

a s = 0.011«;““ (2.26) 

ks = 0.029«;»“, (2.27) 
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within 1.5 a statistical accuracy, which first additionally includes 


CXg > 0 


(2.28) 


possibility alongwith large running compared to the Planck (2013) data. 

At the next to leading order, the simplest way to modify the power law param¬ 
eterization is to incorporate the effects of higher order Logarithmic corrections in 
terms of the presence of non-negligible running, and running of the running of the 
spectral tilt as shown in Eq (2.17), which involves higher order slow-roll corrections 9 . 

After substituting Eq (2.17) in Eq (2.14), I will show that additional information 
can be gained from our analysis: first of all it provides more accurate and improved 
bound on tensor-to-scalar ratio in presence of non-negligible running and running of 
the running of the spectral tilt. In our analysis super-Planckian physics doesn’t play 
any role as the effective theory puts naturally an upper cut-off set by the Planck 
scale. Consequently the prescription only holds good for: 


1. sub-Planckian VEVs, 0 O < M p , 

2. field excursion, A0 < M p for inflation. 

Both of these outcomes open a completely new insight into the particle physics 
motivated models of inflation, which are valid below the Planck scale. 

Further note that the momentum integral has non-monotonous behaviour of the 
slow-roll parameters (ey,rjy) within the interval, k e < k < k cm b, which implies that 
6y and rjy initially increase within an observable window of e-foldings (which I will 
define in the next section, see Eq (3.1)), and then decrease at some point during 
the inflationary epoch when the observable scales had left the Hubble patch, and 
eventually increase again to end inflation [21, 24], 

After substituting Eq (2.17) in the left hand side of I Eq (2.14), I obtain: 



dk r(k ) 

TV “T“ 




a +!K£)+t ln2 T) 


k\ A+Bh <£)+ cln2 (£) 


(2.29) 


where 


Let us substitute, 




c 

12 ' 


k/h = In y, 


(2.30) 


9 It is important to note that when Ref. [25] first derived a bound on large tensor-to-scalar 
ratio for super-Planckian inflationary models (with A0 > M p ), the above mentioned constraints on 
as, ns were not taken into account due to lack of observational constraints. 
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to simplify the mathematical form of the above Eq (2.29). Consequently, 1 get: 


f dk / r (^T fr(k+) f _dy_ ,, xA+B\n{\ny)+Cln 2 {\ny) (2 31) 

J K fcV 8 V8 J eke/ ^ylny [ V) 

To evaluate the integral analytically, 1 apply the following technique. Let us consider: 



(lny)", where a << 1 

(2.32) 

where the exponent a is 

defined as: 



a — A + B In (In y) + C In 2 (In y) 

(2.33) 

where 

1^1, \B\, \C\ « 1 

(2.34) 

with 

\A\>\B\>\C\. 

(2.35) 

Now, for a « 1, which 

is typically the case, one can expand the function mentioned 


in Eq (12) as 10 : 

(lny)" = 1 + a ln(hiy) + • • • (2.36) 


Let us take first two terms in the right hand side of the series expansion. This 
finally results in: 


" k * dk Ir(k ) 
k 


r(h) 


dy 


{1 + [A + D In (In y) + C In 2 (In y)] ln(ln y )} , 


r(h 


r(h 


J e k e /ki y In y 

[(1 — A + 2B — 6 C) In y + (A — 2B + 6C) (lny) ln(lny) 
+ (B- 3 C) (lny) ln 2 (lny) + C'(lny)(ln(lny)) 3 ]^ , 

h 


(1 - A + 2B - 6 C) 




- (A-2B + 6 C) ^ In ( ^ 




— (B — 3C) In 2 (\ -C^ln 3 


r(k* 


„ a b c 

2 — — T — — — 

2 2 2 


h 

ke 

K 


k. 






- --- + r-l win w 


kf 


c\ k, 


--— hT —-- In — 

4 J An \kj 12 An 


c k, 


(2.37) 


10 One can verify the validity of a « 1 for a generic single field slow roll inflation, within the 
interval 8.2 x 1CP 11 Mpc^ 1 < k < 0.056 Mpc -1 . 
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In any arbitrary momentum scale k spectral tilt, running of the tilt, and running 
of the running of the tilt for the scalar and tensor perturbations can be written as: 


ns(k) -1 = 


n T {k) 

as,r(k) 

Ks,T{k) 


din P s (k) 
d In k 


ks(K) 


In 2 [rr- + 


= n s (K) - 1 + a s {K) In J + ^ 

_ din P T {k) ( k\ k t {K 

= ~dwT = nr(K) + aT{h) ,n [kJ + — 


(2.38) 


k 


,n Uj + '< 239) 

dnsrik ) d 2 lnPsT(k ) . ,, . , / k \ 

= a ^> + 111 U + ''' (2 ' 40) 


_ das,r(k) d 2 ns,T(k) d 3 In Ps,r{k) 


din k 


d In k 2 


d In k 3 


Ks,T(k*) H-• 


(2.41) 


Also at scale k tilt, running ang running of the running in tensor-to-scalar ratio 
can be expressed as: 


n r (k ) 
a r (k) 


dr{k ) 


K r (k) = 


= r(k ) 


din k 
dn r (k ) d 2 r(k) 
dink d In /c 2 


a + Mill + ?ln 2 ) +••• 


= [an r (fc) + 6r(/c)] + [bn r (k) + cr{k)\ In — + 


k 


K 


cn r (k) 


(2.42) 


In 2 — + • • ■ (2.43) 


_ da r (k ) d 2 n r (k) d 3 r{k) 


din k 


d In k 2 d In k 3 


[aa> r (k) + 2 bn r (k) + cr(/c)] + [&a r (fc) + 2cn r (/c)] In ( — )+••• 


k 


(2.44) 


Here at an arbitrary scale k, the parameters a, b and c defined as: 


a{k) = 


d\nr{k ) 
din k 


o + Hnld+bn 2 f4| + ... 


da(k ) d 2 lnr(/c) 

o(fc) = 


din k 


dink 2 


b + c In ( — + 


_ da r (k ) d 2 n r (k) d 3 lnr(k ) 

c(k) = 


d In k d In /c 2 


din A: 3 


c+--- 


(2.45) 

(2.46) 

(2.47) 


In the present context the potential dependent slow-roll parameters: (ey, rjy, • • •), 
satisfy the joint Planck (2013)+WMAP-9 constraints, which imply that [8]: 


ey < 10 (within 1.5cr C.L .), 
5 x 10~ 3 < \r]v\ < 0.021 (within 1.5a C.L .), 


(2.48) 

(2.49) 
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and also satisfy the joint Planck (2015) I WMAP-9 I high-L (TT) constraints, which 
imply that [15]: 


e y < 0.011 

(within 2cr C.L .), 

(2.50) 

8 x 10~ 3 < \t)v | < 0.021 

(within 1.5 a C.L), 

(2.51) 

8 x 10~ 3 < \&\ < 0.021 

(within 1.5 a C.L .), 

(2.52) 


for which the inflationary potential is concave in nature for both the cases. In the next 
section, I will discuss model independent bounds on the coefficients (V ( 0 *), 1 /'( 0 *), • • •) 
for a generic sub-Planckian VEV inflationary setup, for which I will satisfy the 
joint constraints from: Planck(2013 & 2015)+WMAP-9+high-L, Planck (2015)+BI- 
CEP2/Keck Array and Planck (2014)+WMAP-9+high L+BICEP2 (dust) data sets. 


3 Constraining the scale of effective field theory inflation via 
field excursion 


The number of e-foldings, N{k), can be expressed as [32]: 


N(k) ss 71.21 - In 



+ 




+ 7 ln 

4 



1 3 Wint i 

H- 7 -- In 

12(1 + Wi nt ) 



(3.1) 

where p e is the energy density at the end of inflation, p r h is an energy scale during 
reheating, A;* = a*77* is the present Hubble scale, V* corresponds to the potential 
energy when the relevant modes left the Hubble patch during inflation corresponding 
to the momentum scale Ac*, and Wi n t characterises the effective equation of state 
parameter between the end of inflation, and the energy scale during reheating. 

Within the momentum interval, k e < k < k±, the corresponding number of 
e-foldings is given by, AN, as: 


AN 





+ 




(3.2) 


where (a*, 77*) and (a e 77 e ) represent the scale factor and the Hubble parameter at 
the pivot scale and end of inflation, and I have used the fact that 


■2 V(4>) 

" 3 Ml 


(3.3) 


which is true within the framework of Einstein’s General Relativity. We can estimate 
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the contribution of the last term of the right hand side by using Eq (7.27) as follows: 


In 



In 


In 


a + 0(0* - 0 O ) + 7(0* - 0 O ) 3 + k :(0* 

a + 0(0 e - 0 O ) + q(0 e - 0o) 3 + «(0e 




V 


3 cc \Mp / 


[ 1 +. 


«i 


0o) 4 + --^ 

0O) 4 + --J’ 



«1 


where 

(A0/M p ) « 1, 

and additionally I assume that 


(/ 3M p /a ) << 1. 

Consequently, Eq (3.2) reduces to the following simplified expression: 


(3.4) 


(3.5) 


(3.6) 


AN m In « In ^ ~ 0(8 — 17) efolds . (3-7) 

Within the observed limit of CMB distortion, i.e. AN m 17, the slow-roll parame¬ 
ters, see Eqs. (10.12, 10.13) of Appendix, show non-monotonic behaviour, where the 
corresponding scalar and tensor amplitude of the power spectrum remains almost 
unchanged 11 . Substituting the results obtained from Eq. (10.12) and Eq. (10.13) 
(see Appendix), and with the help of Eq. (3.2), up to the leading order, 1 obtain: 


n =0 



+ ... 


where for AN = 17 and AN = 81 have used 


(3.8) 


(k e /h) « exp(-AlV) = exp(—17) « 4.13 x 10~ 8 , (3.9) 

(kjK) ~ exp(—AiV) = exp(—8) ~ 3.35 x 10 -4 . (3.10) 

Here I will concentrate on 

a, b, 0 (3.11) 

11 In this paper I fix AW ~ 0(8 — 17) e-foldings as within this interval the combined constraints 
from Planck (2013 & 2015)+WMAP-9+high-L and Planck (2015)+BICEP2/Keck Array are sat¬ 
isfied. Additionally I have also mentioned the results by applying the constraints from Planck 
(2014)+WMAP-9+high L+BICEP2 (dust) data to explicitly show that our prescribed methodol¬ 
ogy also holds good for r(fc*) > 0.12. 
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where additionally 


a » b » c 


(3.12) 


case is satisfied (for the details, see [23, 26, 27, 33]). In Eq (3.8) the series appearing 
in the left side of the above expression is convergent, since the expansion coefficients 
can be expressed as: 



( 


( 


\ 


and I have defined a new dimensionless binomial expansion co-efficient (A m ), defined 
as: 



(Vm = 0,1,2,....), (3.14) 


which is obtained from the binomial series expansion from the leading order results 
of the slow-roll integrals stated in the Appendix 12 . 

Note that the expansion co-efficient A m (Vm) are suppressed by, V(<t> o), which is 
the leading order term in a generic expansion of the inflationary potential as shown 
in Eq (1.6) and Eq (7.27)(see also Eq (10.14) in the Appendix). We can expand the 
left side of Eq. (3.8) in the powers of A <j)/M p , using the additional constraint 


A0 < (0 e - 0o) < M p , 


(3.15) 


and I keep the leading order terms in A <f>/M p . 

To the first order approximation - I can neglect all the higher powers of fc e /fc* ~ 
C>(10 -4 — 10 -8 ) from the left hand side of Eq (3.8), within 0(8 — 17) e-foldings of 
inflation [23, 28, 29]. Consequently, Eq. (3.8) reduces to the following compact form 
for r(k+): 



9 

25 


(3.16) 


provided at the pivot scale, k = k* » k e , in this regime 


\w\ » { e VO livi °Vi ''' } 


(3.17) 


12 In Eq. (3.14), and Eqs. (10.12, 10.13) (see Appendix), C p and T) q are Planck suppressed dimen¬ 


sionful (mass dimension M p ( m+2 ) ) binomial series expansion coefficients, which are expressed in 
terms of the generic model parameters (E(0*), E 7 ^*), • • •) as presented in Eq (7.27). 
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approximation is valid. Our expression, Eq. (3.16), shows that large value of r(k+) 
can be obtained for models of inflation where inflation occurs below the Planck 
cut-off. Once the field excursion, \A(p\/M p , and r]y , are known from any type of 
sub-Planckian inflationary setup, one can easily compute the tensor-to-scalar ratio 
by finding the roots. 

The expression, as stated in Eq (3.16), is in the form of a simple algebraic (cubic) 
equation. In order to find the roots of tensor-to-scalar ratio r in terms of the held 
excursion \A(f>\/M p , one has to solve a cubic equation: 


a : 3 + Ux — W = 0, 


(3.18) 


where 


x : 

U : 
W : 


r(h) 

0.27 ’ 

800 , 

27 {VV + 2) ’ 
40000|A0| 

243 M„ ' 


(3.19) 

(3.20) 

(3.21) 


The three roots X\, £ 2 , £3 are explicitly given by: 


x 1 


X2 


x 3 


(i) 1/3 e , Y_ 

Y ^ 2 1 / 3 3 2 / 3 ’ 

(1 + iV 3) U (1 - iV 3) Y 

2 2/3 3 i/3y 2 4 / 3 3 2 / 3 ’ 

(1 - iV 3) U (1 + iV 3) y 

2 2/3 3 i/3y 2 4 / 3 3 2 / 3 ' 


(3.22) 

(3.23) 

(3.24) 


where the symbol, 

y = (9W + V4f/ 3 + 271E 2 ) 1/3 . (3.25) 

Here the complex roots £2 and £3 are physically redundant. The only acceptable 
root is the real one, i.e. X\ which will finally contribute for the estimation of held 
excursion for a specihed value of tensor-to-scalar ratio r{k+). 

Further note that our formulation will also hold true if inflation were to start at 
the hill-top, such as 

0 = 0. (3.26) 


However in this case one would have to proceed similarly by expanding the potential 
around 


0o = 0, 


(3.27) 


and then follow the algorithm I have provided here. 
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Now, it is also possible to recast a(k), b(k), c(k), in terms of r(k), and the slow 
roll parameters by using the relation, Eq. (10.8) (see Appendix): 


a(k * 


r(k * 


- 2 rj V (h) H- 


b(h) « [l6ey(/c*) - 12e v (k ic )rj v (k ir ) + 2£y(kfl) + • 
c(fc*) ~ [—2cry + ■ • • ] , 




(3.28) 


where “ • • ■" involve the higher order slow roll contributions, which are negligibly 
small in the leading order approximation. 

Now at the pivot scale, k * = 0.002 Mpc~ l , the scale of inflation R 1//4 (0*) can be 
expressed as: 

V l, \^)= Qp S (fc*)r(fc*)) ' x v^M p < (2.40 x 10 16 GeV) x ' • (3-29) 

The equivalent statement can be made in terms of the upper bound on the numerical 
value of the Hubble parameter at the exit of the relevant modes: 


EL < 1.38 x 10 14 x 


r{K) 

0.27 


GeV. 


(3.30) 


Combining Eqs. (3.16), (3.30) and (3.29), 1 can now obtain a closed relationships: 


M 

Mr, 


< 


y/V* 


(1.55 x 10- 2 M p ) 2 


V* 


wih) 


(2.78 x 10 _2 M P ) 4 


(3.31) 


\M H> _ Hi Vv (K) 

M p ~ (1.39 x 10“ 4 M p ) (1.99 x 10~ 7 M 2 ) 2 


(3.32) 


where additionally 

\w\ » T lvi a vi ''' } (3.33) 

are satisfied. Similar expressions were derived in Ref. [23] for inflection-point model 
of inflation, where an additional constraint 


ld"(0 o ) = 0 


(3.34) 


is considered for the computation. The above Eqs. (3.31, 3.32) characterize the 
bounds on A</> for: 


0o < M p , (3.35) 

A0 < M p . (3.36) 
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Our conditions, Eqs. (3.16, 3.31), provide new constraints on model building for 
inflation within particle theory, where the inflaton potential is always constructed 
within an effective field theory with a cut-off. Note that 

\w(K)\ >0 (3.37) 

can provide the largest contribution, in order to satisfy the latest bound on the 
tensor-to-scalar ratio, the shape of the potential has to be concave in nature. Fur¬ 
ther applying this input in Eq (3.29), one can get a preferred bound on the scale of 
sub-Planckian VEV inflation as: 

Planck (2013)-fWMAP-9+high L: 


</K < 1.96 x 10 16 GeV 

for 

r* < 0.12, 

(3.38) 

77* < 9.20 x 10 13 GeV 

for 

r* < 0.12, 

(3.39) 

Planck (2014)+WMAP-9+high L+BICEP2 (dust): 


2.07 x 10 16 GeV < tfv* < 2.40 x 10 16 

GeV 

for 0.15 < r* < 0.27, 

(3.40) 

1.03 x 10 14 GeV < 77* < 1.38 x 10 14 GeV 

for 0.15 < r* < 0.27. 

(3.41) 

Planck (2015)-(-WMAP-9+high L(TT): 




</K < 1.92 x 10 16 GeV 

for 

f * < 0.11, 

(3.42) 

77* < 8.80 x 10 13 GeV 

for 

T * < 0.11, 

(3.43) 

Planck (2015)-1-BICEP2/Keck Array : 




a/K < 1-96 x 10 16 GeV 

for 

f * < 0.12, 

(3.44) 

77* < 9.20 x 10 13 GeV 

for 

r* < 0.12. 

(3.45) 

4 Non-monotonic behaviour of slow-roll parameters within 


effective field theory 


Let us now discuss the non-monotonous features of slow-roll parameters ey,r]v, a v 
appearing in the present context of the paper:- 


• To violate the Lyth bound mainly ey and also the other slow-roll parameters 
r )v,£v,o'v must decrease at some point during the inflationary epoch. This is 
evident from the definition of no. of e-foldings, 


N = 


i dcf) 

m p hcmb 



(4.1) 
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(0)1 «v(0) 


ev 0 vs 0-i 


M) I vs (f 



(Ho) illM p (f-MiiM, 

(c) ICyl vs 0- </>o (d) |<4| vs </>- (j)Q 

Figure 1. Non-monotonous evolution of the slow roll parameters are shown with respect 
to (f> — cj) 0- The upper and lower bounds are set by Eqs. (5.11-5.15). 
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(a) ty vs N 


(b) \ t ) V \ vs N 


|^ 2 y| vs N plot 


|g\| vs N plot 



N 


N 


(c) 41 vs N 


(d) |4| vs N 


Figure 2. Non-monotonous evolution of the slow roll parameters are shown with respect 
to number of e-foldings N. The upper and lower bounds are set by Eqs. (5.11-5.15). 
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where, if ey decreases, AN will increase for the same field excursion, A <fi. 
Unfortunately this alone is not sufficient enough to successfully evade the Lytli 
bound and match all observations. 

• At CMB scale, 

k = k cmb (= a,H ), (4.2) 

at least ey must be large enough to generate an observable value of tensor-to- 
scalar ratio r. In our methodology I have considered that the other slow-roll 
parameters r]y, cry also be sufficiently large to confront observation at CMB 
scale. 

• Also ey and the other slow-roll parameters must increase over the AN ps 17 e- 
fold observational window [23, 28, 29]. This is dictated by a combination of the 
spectral index constraint and the observed value of <jg from large-scale structure 
(LSS) [34], which means the spectrum must decrease over the observational 
window. When running of the spectral index is allowed, the best fit value for 
the spectral index is indeed 

n s > 1, (4.3) 

but the running is 

a s < 0, (4.4) 

therefore ey must increase eventually. Also, the value of cr 8 measured inde¬ 
pendently from LSS strongly favours a primordial spectrum that decreases 
amplitude between Planck and LSS scales. 

• After observable scales have left the horizon (k > aH), ey and the other slow- 
roll parameters must quickly decrease. The quick decrease of ey and as well 
as other slow-roll parameters are necessary to generate enough e-folds required 
for inflation. If instead ey decreases gradually within the present context, then 
it will need to eventually decrease to a much smaller value because in such a 


situation 

ey oc (A(p/2M P AN) 2 , 

(4.5) 

T] V oc 2ey + (1/a/2A N), 

(4.6) 

£y oc 2ey [l + 2ey + (y/2/A A)] + (l/v^AA) 2 , 

(4.7) 

cry oc 2ey |4ey 1 + ey + (V2/AN) + (l/\/2A N) 

•J2 + (l/2f,,AA0 2 ] + sa/v^AJV) 2 } , 


(4.8) 

and I finally require 

A 0 < M p 

(4,9) 
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to violate the Lyth bound in the present context. Consequently, the slope of 
ey, rjv, />, (Ty in fig (1) and fig (2) can be computed as: 


Aey/A0 oc (-/ ey/MpAN ), (4-10) 

Ariv/A <t> oc 2(V^/M P A1V) + (l/2M p ^2i/(A!V) 2 ), (4.11) 

A&/A <f) oc 2{^/M p AN) [l + 4e y + (y/2/AN)] + (^/M p (AiV) 2 ) 

+ (1/4M P ^/(A1V) 3 ), (4.12) 

Acry/A(j) oc (ay/Mp^/e^AN) + {%e v ^/e^/M p AN) 1 + e v + (a/2/ A A) 
TSe 2 , [(y/i//M p AA) + (l/M pV ^/(AA) 2 )] 
+(/f//M pV / 2(AiV) 2 ) [V2 + (l/2e y AA) 2 
+ (l/2M p (2ey) 3/2 (A!V) 4 ) + (6^/4M p (AlV) 3 ), 


and 

|Ae y /AA| oc (4ey/A N), 

\Ar) V /AN\ oc (8ey/AlV) + (l/\/2(AA) 2 ), 
|A£ 2 /A1V| oc (8ey/AA) [l + 2e y + (y/2/AN) 


+ 26 


V 


(4.13) 


(4.14) 

(4.15) 

(8ey/AlV) + (V2/(A A) 2 ) 
+ (1/AiV) 3 , (4.16) 


\A(jy/AN\ oc {2oy/ AN) + ( 32e 2 v /AN ) 1 + ey + (y/2/AN) 


+8 e 2 v (4ey/AA) + (a/2/ (AA) 2 ) 


-(v / 2ey/(AA) 2 ) [v / 2 + (l/2eyAA)' 
-(/2/(AA) 4 ) + (6v / 2ey/(A A) 3 ). 


(4.17) 


• Just before the end of inflation ey and the other slow-roll parameters must 
eventually increase again. 

• The change in ey and other slow-roll parameters must also not be too sharp in 
order not to violate slow-roll condition within the prescribed setup. 

• In case of usual monotonic case the end of inflation (</> e ) is fixed by the condition: 

max^ e [ey, |r/y|, |£y|, |n/|] = 1. (4.18) 

But in the present case the non-monotonicity of the slow-roll parameters cannot 
be able to push the maximum field value upto unity at the end of inflation. In 
this case fig (1) suggest that 

max^^ [\vv\] > max 0=< ^ [ey] = max^^ [|cr/|] > max^ e [|£y|] << 1 

(4.19) 
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and this also implies that from the maximization of rjy one can compute 
the field value corresponding to the end of inflation. Eq (4.19) can also be 
represented in the from the following constraint conditions: 


\V"(<Pe) | 

« 

V'fa) 

« 

\V\<t>e)\ 

« 


= 

\V”"(<t>e)\ 

» 

\V\^e)\ 

» 

\V(d>e)\ 

» 


V(<f> e ) 

Ml ’ 

V2V(</>e) 

Mp ’ 

v/2|E(0 e )E"(0 e )|, 

2 My" fa), 

2 Ml 

2M p 2 |E'"(0 e )|, 

M 2 P \V"fa)\. 


(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 


Further substituting Eq (1.6) in Eq (4.23) finally one can End the field value 
(</> e — </>())■ 


The most demanding aspect of the conditions is that the single scalar poten¬ 
tial in a non-trivial fashion reduces ey, r]y, £y, cry after observable scales leaving the 
horizon at k > aH and subsequently in a non-standard way increases ey , r ] V , , ay 
just before the end of inflation. If inflation is, instead, brought to an end by another 
non-trivial mechanism, such as a hybrid transition [22], then one can relax the fifth 
condition in the present context. This makes the task of reconstructing the infla¬ 
tionary potential simpler. However, in this work I have utilized almost all of the 
conditions to reconstruct the structure of inflationary potential. 

Most importantly, the condition that ey first increases and then decreases is 
not very uncommon in the context of inflationary model building and can easily be 
achieved by simply adding a constant vacuum energy correction term to a potential 
that already supports inflation and has increased ey as well. Before the constant 
vacuum energy correction term becomes completely dominant, ey and the slow- 
roll parameters will continue to increase as before. Then, when the constant term 
does come to completely dominate, ey as well as the other slow-roll parameters 
approach towards zero. If the contribution from such vacuum energy dominated 
correction term is sufficiently large enough, this transition will occur before the end of 
inflation. This also suggests that to successfully evade the Lyth bound in the present 
context, the occurrence of this special feature must coincide with the end of the 
observational window and occur with the right magnitude, which is not guaranteed. 
However, the non-monotonic variation of ey as well as the other slow-roll parameters 
are the examples of this physical mechanism, which will take care of all these issues 
successfully. 
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Additionally it is important to note that with non-monotonic evolution of ey will 
modify the power-law feature in the primordial density perturbations. Apart from 
an observable tensor-to-scalar ratio, such non-montonicity criteria of the slow-roll 
parameters generically exhibit the following unconventional features in the present 
context: 


• Existence of a non-negligible, scale-dependent running (ag) and running of the 
running (/eg), of the scalar spectral index (ns). 


• Appearance of a significant increase in the primordial power spectrum on very 
small momentum scales. 


The non-negligible running and running of the running arises because there must 
be significant evolution of ey as well as the other slow-roll parameters while observable 
momentum scales are crossing the horizon at k — aH . This is possible in principle 
for all the evolution of ey to occur only after the observable scales have crossed 
the horizon, it requires sharp features in ey and after that the decrement in ey 
and other slow-roll parameters begin, the greater the decrement must be in order 
to get sufficient e-folds to achieve inflation. If this decrement is huge and very 
faster, then the perturbative approach of the scalar fluctuations, and consequently 
the consistency of inflation, breaks down in the context of inflationary reconstruction 
of potential. This problem can only be addressed by a concurrent sharp fall in 
the energy scale associated with inflation. For details see ref. [35] where such a 
sharp change would be described as an inflationary model with multiple periods of 
inflation, rather than one single period with an evolving ey. Also to explain the non¬ 
monotonicity of slow-roll parameters the reconstruction technique is implemented in 
order to impose a sufficient running and running of the running at the pivot scale. 
The necessary evolution of ey required for consistency of inflation pointing towards 
the fact that the running and running of the running will be scale-dependent. This 
implies that the description of the scalar primordial power spectrum using the simple 
power-law parameterization in momentum scale is not sufficient enough to explain 
non-monotonic evolution of ey and as well as other slow-roll parameters to evade the 
Lyth bound. Finally, the increase in the primordial power spectrum at very small 
momentum scales is an outcome of the necessary decrement of ey and other slow-roll 
parameters outside the observational window. If the amplitude of the primordial 
power spectrum on these scales is sufficiently large enough then primordial black 
holes (PBHs) will be formed 13 . 


13 Present cosmological constraints on PBH’s restrict P$ < 10 2 


See refs. [36-38] for details. 
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5 Reconstruction technique of the structure of inflationary 
effective potential 

Let us now mention the crucial steps followed during the reconstruction of inflationary 
potential:- 

• STEP I : 

First of all 1 need the required data for inflationary observables from an observa¬ 
tional probe at CMB scale. In our case I use Planck (2013)+WMAP-9+high-L, 
Planck (2014)+WMAP-9+high L+BICEP2 (dust), Planck (2015) +WMAP- 
9+high-L (TT) and Planck (2015)+BICEP2/Keck Array data sets for the nu¬ 
merical estimations. Most importantly, I need to implement the reconstruction 
technique in such a way that I can check the validity and convergence of the 
proposed methodology using any observational probe. 

• STEP II : 

Next I assume the slow-roll paradigm within our prescription. Also I demand 
from the held theoretic requirement that the reconstructed potential is renor- 
malizable and for this reason the Taylor series expansion of the potential is 
truncated at the fourth order derivative term. 

• STEP III : 

Further I construct the scale and the Taylor series expansion co-efficients of 
inhationary potential from the derivatives upto fourth order at the CMB scale. 

• STEP IV : 

Hence I use matrix inversion technique to determine the Taylor series expansion 
co-efficients of the potential around the VEV </> 0 of the inhationary potential. 
This is only possible when the system (square) matrix computed from the 
difference in the held value of the potential at the CMB scale and at VEV is 
non-singular. 

• STEP V : 

Finally, to determine the various cosmological parameters and to check the 
validity of the reconstruction technique I ht the reconstructed potential with 
the observed CMB angular power spectra from TT, TE, EE and BB mode 
obtained from WMAP9, Planck (2013, 2014, 2015) alog with BICEP2/Keck 
Array joint data. 

Let us now discuss observational constraints on the Taylor expansion coefficients 
(V(4> o), V' (c/)o), • • •) as appearing in Eq (1.6). 
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Let us first write down V(0*), V '(</>*), E"(0*), •• • in terms of the inflationary 
observables (see Appendix, Eqs. (10.4, 10.6, 10.8)): 


= -P s (h)r(k*)ir z Mp, 


3 


V'ifa) = ^P s (kMk,W\l r —M^ 


¥"(</>*) = ^Ps(h)r(h)n 2 \^n s (h) - 1 + v 8 "* y ) Mp, 


3 r(k*) \ 


V"\4>.) = -Ps(kMk*y 
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+ i 2 (dMV +r(M („ s(M -i + !^)) 2 
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x 


3r(/c*)\ 
“8/ 
1 / r(fc*)V 


■ 


CCs(fc*)- 


r(fc*) 


^(n s (W- 1 + ?hh))-6'' r(MV 



(5.1) 


Therefore, for any chosen sub-Planckian VEV of 0*, I can obtain a matrix equa¬ 
tion characterizing the coefficients: V(<f> o), E'(0o), V"(4>o), ■ ■ ■: 


^ i , 9 ^ ^ 

* 2 6 24 
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= 

V"(0*) 
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^""(0o) , 


v v"'U) 
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(5.2) 


where I have dehne 


$* := (0* - 0o) < M p . 
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The square matrix marked by the symbol is nonsingular, since its determinant is 
nonzero, for which the matrix inversion technique is applicable in the present context. 
Finally, I get the following physical solution to the problem: 
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(5.3) 


For a model independent constraint on the shape of the potential, the parameter, 

d* ~ 0(W~'M p ), (5.4) 


A0~ O(10 _1 M p ), (5.5) 

which is applicable to a large class of sub-Planckian inflationary models to generate 
sufficient amount of inflation within 50 < N to tai < 70. This criteria holds good for 
high scale inflation, but within the regime of sub-Planckian cut-off. 

In order to satisfy the preferred bounds, see Eq. (1.17)-Eq. (1.36), the following 


model independent theoretical constraints on 1/(0*), 1/ (</>*), ■ • ■ have to be imposed: 

Planck (2013)+WMAP-9+high L: 

17(0*) < 0(3.79 - 3.99) x 10“ 9 M 4 , (5.6) 

P'(0*) < 0(4.65 - 4.89) x 10” 10 M 3 , (5.7) 

V" (0*) < 0((—0.41) - 2.42) x 10 _11 Mp, (5.8) 

I/"'(0*) < 0(8.52 - 35.1) x 10 _11 M P , (5.9) 

17""(0*) < 0(0.39 - 4.76) x 10- 9 , (5.10) 

Planck (2014)+WMAP-9+high L+BICEP2 (dust): 

5.27 x 10~ 9 M 4 < 1/(0*) < 9.52 x 10 ~ 9 M*, (5.11) 

2.45 x 10 _lo M 3 < v' (0*) < 1.75 x 10“ 9 M 3 , (5.12) 

2.41 x 10 _11 M 2 < P"(0*) < 3.25 x 10~ 10 M p 2 , (5.13) 

6.35 x 10 _10 M P < v"' (0*) < 7.56 x 10 ~ 10 M P , (5.14) 

5.56 x 10~ 10 < P""(0*) < 4.82 x 10~ 9 , (5.15) 
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Planck (2015)+WMAP-9+high L(TT): 


P(0*) < 0(3.45 - 3.71) x 10- 9 M p 4 , (5.16) 

P'(</>*) < 0(4.04 - 4.35) x 10 -1 °Mp, (5.17) 

P" (</>*) < 0((—2.93) - 1.08) x 10“ n M 2 , (5.18) 

P"'(0*) < 0((—0.48) - (-3.88)) x 10 _10 M P , (5.19) 

P""(<£*) < 0(5.52 - 5.76) x 10“ 9 , (5.20) 

Planck (2015)+BICEP2/Keck Array: 

P(0*) < 0(3.80 - 3.99) x 10~ 9 M 4 , (5.21) 

P # (0*) < 0(4.65 - 4.89) x 10“ 10 Mp, (5.22) 

P" (</>*) < 0((—0.41) - 2.42) x 10 _11 M 2 , (5.23) 

v'" (0*) < 0(8.52 - 35.1) x 10~ 10 M P , (5.24) 

V"" (0*) < 0(0.39 - 4.76) x 10“ 9 , (5.25) 

Now, substituting the above expressions in Eq (5.3), I obtain model independent 
constraints on the Taylor expansion co-efficients at (ft = (po be. V(4>q), V'(<j) o), • • • as: 

Planck (2013)+WMAP-9+high L: 

V{(j ) 0 ) < 0(3.79 - 3.94) x 10~ 9 M 4 , (5.26) 

V\<f> Q ) < 0(4.66 - 4.88) x 10" 10 M p 3 , (5.27) 

V" (0o) < 0((—1.07) - 1.29) x 10 -11 Mp, (5.28) 

P"'(0o) < 0(8.13 - (-1.25)) x 10 _10 M P , (5.29) 

V'"\<f ) 0 ) < 0(0.39 - 4.76) x 10" 9 , (5.30) 

Planck (2014)+WMAP-9+high L+BICEP2 (dust): 

5.26 x 10~ 9 M 4 < P(0 O ) < 9.50 x 10“ 9 M p 4 , (5.31) 

2.44 x 10~ 10 M p < 1/'(0 O ) < 1.74 x 10' 9 M p 3 , (5.32) 

2.10 x 10 _11 M p < V"{(j > 0 ) < 3.22 x 10" 10 A/ p 2 , (5.33) 

6.29 x 10~ 10 M P < V" {(f > 0 ) < 7.08 x 10“ 10 M P , (5.34) 

5.56 x 10~ 10 < V""{(f > 0 ) < 4.82 x 10~ 9 , (5.35) 

Planck (2015)+WMAP-9+high L(TT): 

V{<f > 0 ) < 0(3.41 - 3.67) x 10“ 9 M 4 , (5.36) 

V\(f)o) < 0(4.06 - 4.31) x 10 _1 °Mp, (5.37) 

v" (0o) < 0(0.31 - 7.84) x nr n M 2 , (5.38) 

v"'{<f> o) < 0((—6.00) - (-9.64)) x 10 _1 °Mp, (5.39) 

V""{<f> 0 ) < 0(5.52 - 5.76) x 10~ 9 , (5.40) 
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Planck (2015)+BICEP2/Keck Array: 


V{(p 0 ) < 0(3.75 - 3.95) x 10“ 9 M 4 , (5.41) 

V\(t ) 0 ) < 0(4.70 - 5.03) x 10” 10 M 3 , (5.42) 

V"(<f> 0 ) < 0((- 8.74) - 30.30) x 10' u M p 2 , (5.43) 

v"' (0o) < 0(8.13 - 30.34) x 10 ~ 10 M p , (5.44) 

V""(<t > 0 ) < 0(0.39 - 4.76) x 10“ 9 , (5.45) 


Consequently, the slow-roll parameters (ey, rjv, £,y, cry) are constrained by: 

Planck (2013)+WMAP-9+high L: 


£y < 0(7.66 - 7.72) x 10" 3 , 

\Vv\ < 0(6.14 x 10" 3 -0.019), 

|$| < 0(2.34 x 10“ 6 - 0.027), 

| tTy | < 0(1.58 x 10“ 3 -0.019). 

Planck (2014)+WMAP-9+high L+BICEP2 (dust): 

(5.46) 

(5.47) 

(5.48) 

(5.49) 

€y ~ 0(0.10-1.69) x 10" 2 , 

(5.50) 

| i]y | ~ 0(4.57 x 10" 3 - 0.030), 

(5.51) 

\&\ ~ 0(5.60 x 10" 3 -0.014), 

(5.52) 

| <Ty | ~ 0(2.28 x 10~ 4 -0.017). 

(5.53) 

Planck (2015)+WMAP-9+high L(TT): 


e y < 0(7.02 - 7.03) x 10" 3 , 

(5.54) 

\t) V \ < 0(2.94-8.59) x 10“ 3 , 

(5.55) 

|£ 2 |< 0(1.52 xl0~ 3 -0.012), 

(5.56) 

Wv\ < 0(0.022 - 0.023). 

(5.57) 

Planck (2015)+BICEP2/Keck Array: 


e y < 0(7.70 - 9.67) x 10" 3 , 

(5.58) 

\riv\ < 0(5.20 x 10“ 4 -0.160), 

(5.59) 

\e v \ < 0(5.60 xlO- 3 - 0.123), 

(5.60) 

\cry\ < 0(1.60 x 10“ 3 - 0.023). 

(5.61) 


Further, by applying the joint constraints from Planck (2013) I WMAP-9 | high L, 
Planck (2014) | WMAP-9 I high L+BICEP2 (dust), Planck (2015)+WMAP-9+high L(TT) 
and Planck (2015)+BICEP2/Keck Array I obtain the following model independent 
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bound on field excursion \A(j)\/M p by using Eq (3.16) or Eq (3.31): 

Planck (2013)+WMAP-9+high L: 


< 0(0.239 - 0.241) for r* < 0.12. (5.62) 

Planck (2014)+WMAP-9+high L+BICEP2 (dust): 

0.242 < < 0.354 for 0.15 < r* < 0.27. (5-63) 

~ M p ~ “ 

Planck (2015)+WMAP-9+high L(TT): 

< 0(0.230 - 0.231) for r* < 0.11. (5-64) 

M p ~ v ' 

Planck (2015)+BICEP2/Keck Array: 

< 0(0.223 - 0.242) for r* < 0.12. (5-65) 

Finally following the present analysis I get the following constraints on the tensor 
spectral tilt, tit, running of the tensor spectral tilt, ar, the running of the tensor- 
to-scalar ratio n r and running of the running of tensor spectral tilt Kt and tensor- 
to-scalar ratio tz r as: 

Planck (2013)+WMAP-9+high L: 


n T < 0((—0.0151) - (-0.0154)), (5.66) 

a T < 0((—2.69) - (-4.815)) x 10" 4 , (5.67) 

\n r \ < 0(3.00 - 3.23) x 10~ 4 , (5.68) 

k t < 0((—5.92) - (-32.20)) x 10“ 5 , (5.69) 

K r < 0((—6.86) - 14.53) x 10" 4 , (5.70) 

Planck (2014)+WMAP-9+high L+BICEP2 (dust): 

- 0.019 < n T < -0.033, (5.71) 

-2.97 x 10" 4 <a T < 2.86 x 10" 5 , (5.72) 

2.28 x 10“ 4 < \n r \ < 0.010, (5.73) 

-0.11 x 10“ 4 <k t < -3.58 x 10“ 4 , (5.74) 

-5.25 x 10~ 3 < n r < -6.27 x 10~ 3 , (5.75) 
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Planck (2015)+WMAP-9+high L(TT): 


n T < C>((—0.0140) - (-0.0142)), 

(5.76) 

a T < 0((- 2.98) - (-5.09)) x 10" 4 , 

(5.77) 

\n r \ < 0(2.47 x 10“ 3 ), 

(5.78) 

k t < £>((- 0.07) - 3.46) x 10” 4 , 

(5.79) 

K r < 0((- 0.14) - 2.89) x 10" 3 , 

(5.80) 

Planck (2015)-t-BICEP2/Keck Array: 


n T < £>((—0.0153) - (-0.0154)), 

(5.81) 

a T < £>((—2.69) - (-4.82)) x 10" 4 , 

(5.82) 

\n r \ < 0(2.11-3.75) x 10“ 3 , 

(5.83) 

k t < 0((—1.79) - (-4.72)) x 10" 4 , 

(5.84) 

K r < 0((—0.20) - (-4.49)) x 10“ 3 , 

(5.85) 


Further, if I had set i? to a slightly larger value, d ~ 0(lO^ 1 M p ) as A cj) ~ 
0( 10~ 4 M p ), then the order of magnitude of the numerics would not change, but 
the numerical pre-factors would slightly change. 


6 Higher order consistency relationships in effective theory 


Let us now provide the new set of consistency relationships between slow roll param¬ 
eters for sub-Planckian models of inflation: 
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(6.3) 
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32 



One can compare these relationships with respect to super-Planckian models of infla¬ 
tion where the slow roll parameters vary monotonically, see Re. [39]. Observationally, 
now one can differentiate sub versus super Planckian excursion models of inflation 
with the help of the above consistency relationships. In particular, the slope of 
the tensor modes, see Eq. (6.1), will play a crucial role in deciding the fate of the 
sub-primordial inflation in the early universe. 

In Fig. (1) and Fig. (2), I have shown the evolution of ey, \vv\, |Zq 2 |, Wv\ (see 
Eqs. (10.3)) with respect to \<p — 0 O I and number of e-foldings N respectively, where 
the upper and lower bounds are given by Eqs. (5.26-5.45). In particular, note that 
the evolution of ey is non-monotonic for sub-Planckian inflation, which is in stark 
contrast with the Lyth-bound for the super-Planckian models of inflation, where ey 
can evolve monotonically for polynomial potentials [39]. In future the data would 
be would be sufficiently good to compare the running of the gravitational tensor 
perturbations, nr, for sub-vs-super-Planckian excursions of the inflaton. 

In Fig. (3), Fig. (4(a)-7(c)) and Fig. (8(a)-ll(d)), I have shown the variation of 
In N(k), Ps(k), ns(k), as(k), Pr(k) and r(k ) (applying all the previously mentioned 
joint constraints), at any arbitrary momentum scale k. Here the black dotted line 
corresponds to 

kmax = 0.3 Mpc" 1 for l max = 2500, (6.7) 

the blue dotted line corresponds to 

k min = 4.488 x 10" 5 Mpc -1 for l min = 2, (6.8) 

and in all the plots violet dashed dotted line represents the pivot scale of momentum 
at 

K = 0.002 Mpc" 1 for Z* ~ 80 (6.9) 

at which 

P s (h) = 2.2 x 10" 9 , (6.10) 

n s = 0.96, (6.11) 

a s = - 0 . 02 , ( 6 . 12 ) 

N(k+) = 63.26. (6.13) 

Within 2 < l < 2500 the value of the required momentum scale is determined by the 
relation [31], 

k„, d ~ (6.14) 

Von 

where the conformal time at the present epoch is given by: 

t]q ~ 14000 Mpc. (6.15) 
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In Fig. (12(a)) and Fig. (12(b)), I have shown the total number of e-foldings 

N, with respect to the field (</> — <fi 0 ) and the field evolution A</> with respect 
to AN respectively. In Fig. (12(a)) the allowed 2 a region is shown as obtained 
from Planck (2013)+WMAP-9+high L, Planck+WMAP-9-fhigh L+BICEP2 (dust), 
Planck (2015) I WMAP-9 | high L(TT) and Planck (2015) +BICEP2/Keck Array 
joint constraints. This also shows that the observational scanning region, AN ~ 17 
is consistent with the bound on the sub-Planckian value of the field excursion as 
obtained in earlier section of the paper. Fig. (12(b)) depicts that at AN ~ 8 and 
AN ~ 17 the numerical value of the field evolution turn out to be A cj) ~ 0.14 M p 
and A(f) ~ 0.34 M p respectively. 

In Fig. (13(a)) and Fig. (13(b)), I have shown the variation of Pg, ns and r 0 .oo 2 , 
at the pivot scale fc* = 0.002 Mpc~ x . The overlapping red patch shows the allowed 
region by Planck (2013)-bWMAP-9+high L, PlanckTWMAP-9+high L+BICEP2 
(dust), Planck (2015)TWMAP-9+high L(TT) and Planck (2015)TBICEP2/Keck 
Array joint constraints. The upper green and lower bounds are set by Eqs. (5.26- 

5.45). 

In Fig (14), I have shown r vs ns at the pivot scale: k+ ~ 0.002 Mpc^ 1 . The 
allowed region is shown by the shaded violet colour, for 0.15 < r* < 0.27 and 

O. 952 < ns < 0.967. The green and yellow lines are drawn for lower and upper bound 
of the constraints derived in Eq (5.46-5.61). We have used the relation between ns 
and r* as mentioned in Eq (10.6, 10.8) in the appendix. 

In Fig. (15(a)-15(c)), I have depicted running of the tensor-to-scalar ratio: n r = 
dr/dink, running of the running of the tensor-to-scalar ratio: K r = d?r/d\nk, 
running of the tensor spectral tilt: = dnr/dlnk, running of the running of 

tensor spectral tilt: Kt = d 2 nx/d\nk vs scalar spectral tilt ng. Shaded violet 
colour region is the allowed overlapping region for Planck (2013)+WMAP-9+high L, 
Planck+WMAP-9+high L-I-BICEP2 (dust), Planck (2015)+WMAP-9+high L(TT) 
and Planck (2015) +BICEP2/Keck Array joint constraints which further constrain 
n r , hi r , otT and Kt within the specified range mentioned in Eq (5.67) and Eq (5.70). 
The green and yellow lines are drawn for lower and upper bound on the constraints 
derived in Eq (5.46-5.51), We have used the relation between and ng as 

mentioned in Eq (6.2) and Eq (6.4). 


7 Example of Inflection point inflation within effective theory 


Now I impose, 


V"{<h) = o, 


(7.1) 
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4.4 


lnN(k) vs Ink plot 


lnN(k) vs Ink plot 



(a) Pt vs k 

lnN(k) vs Ink plot 



(c) r(k) vs k 



(b) Pt vs k 

lnN(k) vs Ink plot 



(d) r(k) vs k 


Figure 3. Here I show the variation of the total number of e-foldings N(k), with respect 
to the momentum scale k. The black dotted line corresponds to k max = 0.3 Mpc -1 for 
Imax = 2500, the blue dotted line corresponds to k m i n = 4.488 X 1CF 5 Mpc -1 for l m ,in = 2, 
and in all the plots violet dashed dotted line represents the pivot scale of momentum at 
k * = 0.002 Mpc -1 for Z* ~ 80, at which for 3(a) Ps(k±) = 2.2 x 1CF 9 , ns = 0.9600, as = 
—0.013 and N(k±) = 63.26 (Planck (2013)+WMAP-9+high-L), 3(b) Ps(fc*) = 2.2 x 10“ 9 , 
n s = 0.9600, a s = -0.022 and JV(fc*) = 63.26 (Planck (2014)+WMAP-9+high-L+BICEP2 
(dust)), 3(c) Pg(fc*) = 2.2 x 10 -9 , ns = 0.9569, as = 0.011 and N(k+) = 63.26 (Planck 
(2015)+WMAP-9+high-L(TT)), 3(d) P S (fc*) = 2.2 x 10~ 9 , n s = 0.9600, a s = -0.022 and 
N(k+) = 63.26 (Planck (2015)+BICEP2/Keck Array). Within 2 < l < 2500 the value of 
the required momentum scale is determined by the relation, k req d ~ [31], where the 

conformal time at the present epoch is r/o ~ 14000 Mpc. 
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P$( k) vs k plot 



(a) P s vs k 


n s (k) vs k plot ff S (k) vs k plot 



(b) ns vs k (c) as vs k 


Figure 4. In 4(a), I show the scalar Power spectrum P s (k) 4(b), I show the scalar spectral 
index n s (k), and in 4(c), I show the running of the scalar spectral index a s (k), with respect 
to the momentum scale k. The black dotted line corresponds to k max = 0.3 Mpc -1 for 
Imax = 2500, the blue dotted line corresponds to k m .i n = 4.488 x 10 -5 Mpc -1 for l m in = 2, 
and in all the plots violet dashed dotted line represents the pivot scale of momentum 
at A;* = 0.002 Mpc -1 for Z* ~ 80 at which Planck (2013)+WMAP-9+high L constraints 
Ps{k+) = 2.2 x 10- 9 , ns = 0.9600, as = —0.02 and N(k*) = 63.26 are satisfied. Within 
2 < l < 2500 the value of the required momentum scale is determined by the relation, 
kreqd ~ [31], where the conformal time at the present epoch is tjq ~ 14000 Mpc. 
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P$( k) vs k plot 



(a) P s vs k 


n s (k) vs k plot ffj(k) vs k plot 



(b) ns vs k (c) as vs k 


Figure 5. In 5(a), I show the scalar Power spectrum P s (k) 5(b), I show the scalar spectral 
index n s (k), and in 5(c), I show the running of the scalar spectral index a s (k), with respect 
to the momentum scale k. The black dotted line corresponds to k max = 0.3 Mpc -1 for 
Imax = 2500, the blue dotted line corresponds to k m i n = 4.488 x 10 -5 Mpc -1 for l rn in = 2, 
and in all the plots violet dashed dotted line represents the pivot scale of momentum at 
= 0.002 Mpc -1 for Z* ~ 80 at which Planck (2014)+WMAP-9+high-L+BICEP2 (dust) 
Ps(k*) = 2.2 x 10- 9 , ns = 0.9600, as = —0.02 and N(k*) = 63.26 are satisfied. Within 
2 < l < 2500 the value of the required momentum scale is determined by the relation, 
kreqd ~ [31], where the conformal time at the present epoch is tjq ~ 14000 Mpc. 
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P$( k) vs k plot 



(a) Ps vs k 

«s( k) vs k plot ff S (k) vs k plot 



Figure 6. In 6(a), I show the scalar Power spectrum P s (k) 6(b), I show the scalar spectral 
index n s (k), and in 6(c), I show the running of the scalar spectral index a s (k), with respect 
to the momentum scale k. The black dotted line corresponds to k max = 0.3 Mpc -1 for 
Imax = 2500, the blue dotted line corresponds to k m i n = 4.488 x 10 -5 Mpc -1 for l m in = 2, 
and in all the plots violet dashed dotted line represents the pivot scale of momentum at 
A* = 0.002 Mpc -1 for ~ 80 at which Planck (2015)+WMAP-9+high-L(TT) constraints 
Ps{k *) = 2.2 x KT 9 , ns = 0.9569, as = 0.011 and IV(A;*) = 63.26 are satisfied. Within 
2 < l < 2500 the value of the required momentum scale is determined by the relation, 
kreqd ~ [31], where the conformal time at the present epoch is tjq ~ 14000 Mpc. 
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P$( k) vs k plot 



(a) P s vs k 


n s (k) vs k plot ffj(k) vs k plot 



(b) ns vs k (c) as vs k 


Figure 7. In 7(a), I show the scalar Power spectrum P s (k) 7(b), I show the scalar spectral 
index n s (k), and in 7(c), I show the running of the scalar spectral index a s (k), with respect 
to the momentum scale k. The black dotted line corresponds to k max = 0.3 Mpc -1 for 
Imax = 2500, the blue dotted line corresponds to k m .i n = 4.488 x 10~ 5 Mpc -1 for l m in = 2, 
and in all the plots violet dashed dotted line represents the pivot scale of momentum at 
A;* = 0.002 Mpc -1 for Z* ~ 80 at which Planck (2015)+BICEP2/Keck Array joint constraint 
Ps(k+) = 2.2 x Hr 9 , ns = 0.96, as = —0.02 and A r (&*) = 63.26 are satisfied. Within 
2 < l < 2500 the value of the required momentum scale is determined by the relation, 
kreqd ~ [31], where the conformal time at the present epoch is tjq ~ 14000 Mpc. 
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Pt( k) vs k plot 


Pt( k) vs k plot 




10- 5 10- 4 0.001 0.01 0.1 1 


io - 5 io - 4 0.001 0.01 0.1 1 


k(inMpc *) 


k(inMpc *) 


(a) Pt vs k 

r(k) vs k plot 


(b) Pt vs k 

r(k) vs k plot 




(c) r(k) vs k 


(d) r(k) vs k 


Figure 8. In 8(a), I show the Tensor spectrum Px{k) by assuming ro .002 ~ 0.12 8(b), 
I show the Tensor spectrum Px(k) by assuming ro .002 ~ 0.12, 8(c), I show the scale de¬ 
pendence of tensor-to scalar ratio r(k) with k by assuming ro .002 ~ 0.12, and in 8(d), I 
show the scale dependence of tensor-to scalar ratio r(k) with k by assuming ro .002 ~ 0.12. 
The black dotted line corresponds to k max = 0.3 Mpc^ 1 for l ma x = 2500, the blue dotted 
line corresponds to k m i n = 4.488 x 10 -5 MpW 1 for l m i n = 2, and in all the plots violet 
dashed dotted line represents the pivot scale of momentum at /c* = 0.002 Mpc^ 1 for Z* ~ 80 
at which Planck (2013)+WMAP-9+high L constraints Ps(k+) = 2.2 x 10~ 9 , ns = 0.9600, 
as = —0.013 and N(k+) = 63.26 are satisfied. Within 2 < l < 2500 the value of the required 
momentum scale is determined by the relation, k req d ~ [31], where the conformal time 

at the present epoch is t]q ~ 14000 Mpc. 
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Figure 9. In 9(a), I show the Tensor spectrum Pj-(fc) by assuming ro .002 ~ 0.27 9(b), I 
show the Tensor spectrum Pr{k) by assuming ro .002 ~ 0.15, 9(c), I show the scale depen¬ 
dence of tensor-to scalar ratio r(k ) with k by assuming ro .002 ~ 0.27, and in 9(d), I show 
the scale dependence of tensor-to scalar ratio r(k ) with k by assuming ro .002 ~ 0.15. The 
black dotted line corresponds to k max = 0.3 Mpc~* for l ma x = 2500, the blue dotted line 
corresponds to k m i n = 4.488 x 10~ 5 Mpc -1 for l m i n = 2, and in all the plots violet dashed 
dotted line represents the pivot scale of momentum at fc* = 0.002 Mpc -1 for Z* ~ 80 at 
which Planck (2014)+WMAP-9+high L+BICEP2 (dust) constraints Ps(k *) = 2.2 x 10 -9 , 
ns = 0.96, as = —0.022 and N(k+) = 63.26 are satisfied. Within 2 < l < 2500 the value 
of the required momentum scale is determined by the relation, k req d ~ [31], where the 

conformal time at the present epoch is r/o ~ 14000 Mpc. 
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Figure 10. In 10(a), I show the Tensor spectrum Px{k) by assuming ro .002 ~ 0.11 
10(b), I show the Tensor spectrum Pr(k) by assuming ro .002 ~ 0.11, 10(c), I show the scale 
dependence of tensor-to scalar ratio r(k) with k by assuming ro .002 ~ 0.11, and in 10(d), I 
show the scale dependence of tensor-to scalar ratio r{k) with k by assuming ro .002 ~ 0.11. 
The black dotted line corresponds to k max = 0.3 Mpc^ 1 for l ma x = 2500, the blue dotted 
line corresponds to k m i n = 4.488 x 10 -5 MpW 1 for l m i n = 2, and in all the plots violet 
dashed dotted line represents the pivot scale of momentum at A;* = 0.002 Mpc^ 1 for Z* ~ 80 
at which Planck (2015)+WMAP-9+high L constraints Ps(k+) = 2.2 x 10~ 9 , ns = 0.96, 
as = 0.011 and N(k+) = 63.26 are satisfied. Within 2 < l < 2500 the value of the required 
momentum scale is determined by the relation, k req d ~ [31], where the conformal time 

at the present epoch is t]q ~ 14000 Mpc. 
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Figure 11. In 11(a), I show the Tensor spectrum Px{k) by assuming ro .002 ~ 0.12 
11(b), I show the Tensor spectrum Pr(k) by assuming ro .002 ~ 0.12, 11(c), I show the 
scale dependence of tensor-to scalar ratio r(k ) with k by assuming ro .002 ~ 0.12, and 
in 11(d), I show the scale dependence of tensor-to scalar ratio r(k ) with k by assuming 
r 0.002 ~ 0.12. The black dotted line corresponds to k max = 0.3 Mpc -1 for Imax = 2500, the 
blue dotted line corresponds to k m i n = 4.488 x 10 - ' 5 Mpc” 1 for l m i n = 2, and in all the plots 
violet dashed dotted line represents the pivot scale of momentum at &:* = 0.002 Mpc -1 for 
l * ~ 80 at which Planck (2015)+BICEP2/Keck Array joint constraints Pg(fc*) = 2.2 x 10 -9 , 
ns = 0.96, as = —0.022 and N(k+) = 63.26 are satisfied. Within 2 < l < 2500 the value 
of the required momentum scale is determined by the relation, k req d ~ [31], where the 

conformal time at the present epoch is r/o ~ 14000 Mpc. 


43 
















N vs (0—0o) plot 



(a) N vs (<j> - (f > 0 ) 



5 10 15 20 


AN 

(b) A<fi vs AN 

Figure 12. We show in 12(a) the total number of e-foldings N. with respect to the held 
(0 — 4>q) and 12(b) the held evolution A(f>, with respect to AN. Within the depicted pa¬ 
rameter space 2a constraints are allowed as obtained from Planck (2013)+WMAP-9+high- 
L, Planck (2014)+WMAP-9+high L+BICEP2 (dust), Planck (2015)+WMAP-9+high- 
L(TT) and Planck (2015)+BICEP2/Keck Array data. This also shows that the observa¬ 
tional scanning region, AN ~ 17 [23, 28, 29] obtained from CMB distortions is consistent 
with the sub-Planckian value of the held excursion. Here | A(p\ = \ (cf> crn b — </>o) — (</> e — </>o) | = 
| (frcmb 0e| ~ | 0e|- 
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Figure 13. We have shown the variation of (a) P$ vs ns, and (b) Ps vs ro.oo 2 j f° r the pivot 
scale k* = 0.002 Mpc -1 . The overlapping red patch shows the 2cr allowed region by the 
joint constraints obtained from Planck (2013)+WMAP-9+high-L, Planck (2014)+WMAP- 
9+high L+BICEP2 (dust), Planck (2015)+WMAP-9+high-L(TT) and Planck (2015)+BI- 
CEP2/Keck Array data. The upper (green) and lower ( ) bounds are set by Eqs. (5.26- 

5.45). 
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Figure 14. We show the joint ler and 2 a CL. contours using Planck 2013, Planck 2015 
+TT+low P and Planck 2015 +TT+TE+EE+low P data for r vs ns plot at the momentum 
pivot fc* ~ 0.002 Mpc -1 . The small circle on the left corresponds to N = 50, while the right 
big circle corresponds to N = 60. The allowed regions are shown by the shaded violet colour 
for r* < 0.12 and 0.952 < ns < 0.967. The blue and green bubbles and stars are drawn 
for saddle point MSSM inflation and inflection point MSSM inflation which are compatible 
with the constraints derived in Eq (5.46-5.49) and Eq (5.54-5.61). For Eq (5.50-5.54) the 
range of r* is outside the present upper bound from Planck 2015 and BICEP2/Keck Array 
joint data sets. The vertical black coloured lines are drawn to show the bounded regions of 
the sub-Planckian inflationary model, along which the number of e-foldings are fixed. 


in order to study the inflection point scenario. The potential is given by 14 : 

V'"(d > ! V""(ch 1 

V(fl = V(<h) + H0o)(0 - <M + —^(0 - 0o) 3 + —^(0 - 0o) 4 + • • • , (7.2) 


14 The inflection point inflation has been studied in Refs. [40], with a constant potential energy 
density V((f>o). 
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Figure 15. We show the (a) tensor spectral tilt nr, (b) running of the tensor spectral tilt 
a-r = driT/dink, (c) running of the running of tensor spectral tilt kt = d 2 nx / d In k vs scalar 
spectral tilt ns plot. The small circle on the left corresponds to N = 50, while the right 
big circle corresponds to IV = 70. Shaded violet coloured regions are the allowed regions 
for Planck (2013)+WMAP-9 high L, Planck - WMAP-9+high L+BICEP2 (dust), Planck 
(2015)+WMAP-9+high L(TT) and Planck (2015) -fBICEP2/Keck Array joint data sets 
which will further constrain ax and nx within the specified ranges mentioned in Eq (5.67) 
and Eq (5.70). The green and yellow lines are drawn for lower and upper bounds on the 
constraints derived in Eq (5.46-5.61), The vertical black coloured lines are drawn to show 
the bounded regions of the sub-Planckian inflationary model, along which the number of 
e-foldings are fixed. 
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We can express U(</>*), W(0*), • • • in terms of 0 O for a sub-Planckian regime as: 

■) 9 3 ) 9 4 

m) = fo(0 o )+ *y{<h) + f fo"'(0 o ) + gfo""(4), 

,92 „q 3 

= fo'(0o) + f fo'"(0o) + xffo""(0o), 

2 6 

fo" (</>*) = *y"{!fo) + ffo""(0o), 

V"'(&) = fo"'(0o)+^fo""(0o), 

V""fa) = fo""(0o). (7.3) 

Using Eq (7.3), I obtain a particular solution for the coefficients: U(0o), V{^> 0 ), • • ■, 
which can be written as 15 : 

i 9 3 

U(0o) = fo(0*) - 0* V(<M + yfo"(</>*) - 

,92 ,q 3 

Vifh) = Vfa) - yU"U) + yfo""(<M, 

U'"(0o) = U"'(0*)-^U""(0*), 

fo""(0o) = fo""(0*). (7.4) 

Now using the bound on U(0*), U'(0*), • • • as mentioned in Eqs. (5.26-5.45), I can 
obtain the following constraints on the coefficients of U(0 O ), V((f) 0 ), ■ ■ •: 

Planck (2013)+WMAP-9+high L: 

U(0 0 ) < 0(3.79 - 3.94) x HU 9 M 4 , (7.5) 

V\<j> 0 ) < 0(4.66 - 4.88) x 10“ 10 M p 3 , (7.6) 

V"'{<j> 0 ) < 0(8.13 - (-1.25)) x 10~ 10 M p , (7.7) 

V"'\(po) < 0(0.39 - 4.76) x 10“ 9 , (7.8) 

Planck (2014)+WMAP-9+high L+BICEP2 (dust): 

5.26 x 10~ 9 M 4 < V(</> 0 ) < 9.50 x 10 ~ 9 M*, (7.9) 

2.44 x 10~ 10 M 3 < V\<t> 0 ) < 1.74 x 10~ 9 M 3 , (7.10) 

6.29 x 10~ 10 M P < V"\<j>o) < 7.08 x 10” 10 M P , (7.11) 

5.56 x 10~ 10 < 0 ) < 4.82 x 10“ 9 , (7.12) 

Planck (2015)+WMAP-9+high L(TT): 

VM < 0(3.41 - 3.67) x 10“ 9 M p 4 , (7.13) 

fo'(0o) < 0(4.06 - 4.31) x 10“ 10 M 3 , (7.14) 

V"\<f> 0 ) < 0((—6.00) - (-9.64)) x 10 _10 M p , (7.15) 

V""(<j> 0 ) < 0(5.52 - 5.76) x 10“ 9 , (7.16) 

15 There will be in general 2 solutions around an inflection point, here I will provide one of the 
two solutions which is the most interesting one for the general case of study. 
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Planck (2015)+BICEP2/Keck Array: 


V(<j> o) < 0(3.75 - 3.95) x 10 ~ 9 M 4 , (7.17) 

U'(0o) < 0(4.70 - 5.03) x 1(T 10 M 3 , (7.18) 

v'" (0o) < 0(8.13 - 30.34) x 10~ 10 M p , (7.19) 

U""(0o) < 0(0.39 - 4.76) x 10“ 9 , (7.20) 


We can compare these results with that of the constraints mentioned in Eq (5.46- 
5.61) for a generic sub-Planckian inflationary setup for = 0* ~ 0o ~ 10 _ 1 M p for a 
sub-Planckian VEV model of inflation. We End a very nice agreement which testifies 
the power of a model independent reconstruction of the potential. 

Let us give an example of Hubble induced supergravity motivated MSSM in¬ 
flation which is guided by inflection point prescription. For the potential under 
consideration, I have 

Ho = 3 H 2 M 2 p ~ M s 4 » mj|0| 2 , (7.21) 

where ~ (P(TeV) is the soft mass. In this case the contributions from the 
Hubble-induced terms are important compared to the soft SUSY breaking mass, 
TO 0 . The potential, after stabilizing the angular direction of the complex scalar held 
0 = |0|exp[i$], see [31, 41, 42], reduces to a simple form along the real direction, 
which is dominated by a single scale, i.e. H ~ 7/*: 


U(0) 


V 0 + c h H 2 \<P\ 2 


a H H(j) n A 2 |0| 2(n_1) 

nM ™~ 3 + Mp (n-3) 


(7.22) 


where I take A = 1, and,the Hubble-induced mass parameter is ch and the trilinear 
A term is a#. Fortunately for this class of potential given by Eq (7.22), inflection 
point inflation can be accommodated, when 


a 2 H » 8 (n - 1 )c H . 


(7.23) 


This can be characterized by a fine-tuning parameter, 5, which is defined as: 


l h 


8 (n - 1 )c H 


= 1 - 


n — 2 


When |<5| is small 16 , a point of inflection 0 O exists, such that 


with 


00 — 


V" (0 O ) = 0, 


l/n-2 

, ) + o{6 2 ' 

[n — 1 ) y 


(7.24) 

(7.25) 

(7.26) 


16 We will consider a moderate tuning of order S ~ 10 4 between ch and an- 
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Figure 16. We show the variation of the n^6 MSSM inflationary potential with respect 
to the field 0, which satisfies all the obtained constraints from recent Planck 2015 data and 
Planck 2015 +BICEP2/Iveck Array joint constraints. 


For S < 1, I can Taylor-expand the inflaton potential around an inflection point, 
0 = 0o, as [31, 41, 42]: 

V{4>) = a + 0(0 - 0 O ) + 7(</> - 0 O ) 3 + «(0 - 0o) 4 H -, (7.27) 


where the expansion coefficients are now given by: 


a = V^o) = V 0 + 


(n- 2) 2 i (n — 2) 
n(n — 1) 


+ 


S 2 \c H H i <& + 0{P), 


P = V’((i >, 0 ) = 2 


n — 2 


5 2 c„77 2 0 o + O(<5 4 ), 


7 = 


K = 


v"\(j) o) c h H 2 f 2 (n — l)(n — 2) 3 


3! 

V"" (0o) 

4! 

c h H 2 

0g 


00 


4(n — 2)" 


h 2 +0(h 4 ) : 


12(n-2) 


3 (n - 1) (n - 2) (n - 3) (7n 2 - 27 n + 26) 


(7.28) 

(7.29) 

(7.30) 

(7.31) 

5 2 ) + 0{5 a ). 


Note that once I specify Ch and 77*, all the terms in the potential are determined. In 
this regard the potential indeed simplifies a lot to study the cosmological observables. 
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As an concrete example, I considered n — 6 case, where the flatness of the 
superfield <f> is lifted by the non-renormalizable operator. This is appropriate for 
both udd and LLe flat directions. 

We fix A = 0(1) and 8 ~ 10 -4 . In order to satisfy the Planck observational 
constraints on the amplitude of the power spectrum, 2.092 x 10 -9 < Pg < 2.297 x 1CP 9 
(within 2cr), spectral tilt 0.958 < ns < 0.963 (within 2cr), sound speed cs — 1 
(within 2er), and tensor-to-scalar ratio r k < 0.12, I obtain the following constraints 
on our parameters for H in f > m$ ~ 0(TeV), where successful inflation can occur 
via inflection point 

c H ~ 0(10 - 10“ 6 ), for 10' 22 < r* < 0.12 

a H ~ 0(30 - 1(T 3 ), for 1(T 22 < r* < 0.12 

M s ~ 0(9.50 x 10 10 - 1.77 x 10 16 ) GeV , for 10~ 22 <r k < 0.12 . (7.32) 

This analysis justifies that the Hubble correction induced MSSM inflation is one of 
the examples of inflection point inflation where the present reconstruction technique 
holds good perfectly. 

8 Example of Saddle point inflation within effective theory 

Now I impose, 


E'(0 o ) = O, (8.1) 

E'U) = 0, (8.2) 

o) = 0, (8.3) 

in order to study the saddle point scenario. The potential is given by 4 ': 

V""((h 1 

1/(0) = V (0o) + 24~ (0 " 00)4 + • • • « ^ 8 - 4 ) 


We can express E(0*), E , (0*), • • • in terms of 0o for a sub-Planckian regime as: 

v(A) = 

v'(M = 

6 

E'"(0*) = ^E""(0 O ), 

V""(0*) = E""(0 O ). (8.5) 

17 The saddle point inflation has been studied in Refs. [43]. 
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Using Eq (8.5), I obtain a particular solution for the coefficients: U(0 O ), V'(<f>o), • ■ ■, 
which can be written as 18 : 

,q3 c,,q4 

V(j> o) = V(0*) - d*U'(0*) + ^U'"(0*) - U (0*), 


V""(0o) = V""(0*). (8.6) 

Now using the bound on U(0*), U'(0*), • • • as mentioned in Eqs. (5.26-5.45), I can 
obtain the following constraints on the coefficients of U(0 O ), U (0o), • • •: 

Planck (2013)+WMAP-9+high L: 

V (0o) < 0(3.79 - 3.94) x 10“ 9 M 4 , (8.7) 

v"" (0o) < 0(0.39 - 4.76) x KT 9 , (8.8) 

Planck (2014)+WMAP-9-f high L+BICEP2 (dust): 

5.26 x 10~ 9 M 4 < V (0o) < 9.50 x 10“ 9 M 4 , (8.9) 

5.56 x KT 10 < U""(0 O ) < 4.82 x 10“ 9 , (8.10) 

Planck (2015)+WMAP-9+high L(TT): 

V (0o) < 0(3.41 - 3.67) x 10“ 9 M 4 , (8.11) 

U""(0o) < 0(5.52 - 5.76) x 10" 9 , (8.12) 

Planck (2015)+BICEP2/Keck Array: 

V (0o) < 0(3.75 - 3.95) x 10" 9 M 4 , (8.13) 

v"" (0o) < 0(0.39 - 4.76) x 10" 9 , (8.14) 


We can compare these results with that of the constraints mentioned in Eq (5.46- 
5.61) for a generic sub-Planckian inflationary setup for $* = 0* — 0o ~ 10 _1 M p for a 
sub-Planckian VEV model of inflation. We find a very nice agreement which testifies 
the power of a model independent reconstruction of the potential. 

Let us consider a concrete example of soft SUSY breaking induced MSSM in¬ 
flation which is guided by the principle of saddle point inflation. Considering the 
contribution from QQQL , QuQd, QuLe and Trade, the n — 4 indexed SUSY flat 
directions the effective potential within the framework of MSSM can be written as 
[43]: 

1 \A A 2 libl 6 

v&, 0) = V o + -ml |0| 2 + —10| 4 cos (4(9 + 9 A ) + (8.15) 

ls There will be in general 2 solutions around an inflection point, here I will provide one of the 
two solutions which is the most interesting one for the general case of study. 
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where Vo is the vacuum energy dominated term which mimics the role of cosmological 
constant, represents the soft SUSY breaking mass term, the inflaton \<fi\ is the 
radial coordinate of the complex scalar field <f> = |0|e* e and the second term is the 
trilinear A-term which has a periodicity of 27T in 2D along with an extra phase 6 a- 
The radiative correction slightly affects the soft term and the position of the saddle 
point in potential valley. One can tune the vacuum energy term Vo ~ 0. But for the 
generality 1 keep this term. Once this term is switched on in the effective potential, 
the scale of the potential goes up to the GUT scale. This will change the value 
of tensor-to-sclar ratio upto r(/c*) ~ 0.12. But the other inflationary observables 
computed from the model is insensitive to the addition of vacuum energy term Vo- 
After minimizing with respect to the angular coordinate 6 the effective potential 
takes the following form [43]: 


U(0) = U 0 + -ml | 


XA 

4AL 


+ 


M 2 


(8.16) 


using which the position for the saddle point is computed from the model as: 


0 o — 

where 

as: 




4A (D 3 + 3) 


a{i+ L Aj± 1 + Y) - (£>, + 1) (D 3 + 3) 


N 

(8.17) 

D,, D 2 , D 3 are the contribution from one loop radiative corrections appearing 


A 

A 

m 2 , 


Ao 


1 + 2 D 3 In 



Alo 

ml 


1 + 2 D 2 In ( 

1 + 2D, In ( — ) 

Vho / 


l + 2D 3 \n(^) 

V/V 


(8.18) 

(8.19) 

( 8 . 20 ) 


Also the trilenear A term and the radiative correction term D 3 satisfy the following 
constraints: 


A 

D 3 


a/2 (D 3 + 3) G 1 G 2 G 37 TI 0 , 
M P A 0 

4^000 (37 + 60 In (^)) 


(l3+ 12 In 



2mlD 1 M p A _ 2OA o 0q A 

A o ^o0o V KAo ) 


( 8 . 21 ) 

( 8 . 22 ) 
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Here G x, G 2 and G 3 is given by: 


G i 
G 2 

g 3 


(■Pi +1) 

(£>3 + 3) 

(D, +1) 

(1 + f) 

(£>3 + 3) 


(15+llA>)-(3Bi + l) , 

~ (3-Di + 1)^1 + 

(11^3 + 15)- ( 3 +^) 


D 2 

2 



(8.23) 

(8.24) 

(8.25) 


After applying the saddle point technique aroun </>o the effective potential can be 
recast as [43]: 


v((j>) = a x + a 2 (0-0o) 4 , 


(8.26) 


where Ax and A 2 are the Taylor expansion co-efficients defined as: 


Ax = V(<j> 0 ) 


y , m l M P 




311+——— 
' 2 6 


1 + 2 D 1 In ( — 
Mo 


— 2 ( 1 + — — — 

' 2 6 


(8.27) 


2 r 


1 + 2D 2 In ( — 
Mo 


A 2 — 


U""(0 O 

4! 


tTln 


D\ D 3 

24v/6 (f>l V + ~2 ~6~ 


360 


2 v / 6T>x 


—= - 12V6 + 684 D 3 - 50 VQD 2 - 

v/6 (i+f-tr 


2 (- 12^6D 2 ) In 

V Vq J Vmq 


(8.28) 


In order to satisfy the Planck observational constraints on the amplitude of the 
power spectrum, 2.092 x 10 -9 < Ps < 2.297 x 10” 9 (within 2cr), spectral tilt 0.958 < 
ns < 0.963 (within 2a), sound speed cs = 1 (within 2cr), and tensor-to-scalar ratio 
r* < 0.12, I obtain the following constraints on our parameters, where successful 
inflation can occur via saddle point. 


Ax ~ e>(10~ 36 - 10“ 9 )M 4 , for 10- 29 < r* < 0.12 

A 2 ~ C>(10' 13 - 10“ 9 ), for 10“ 29 < r* < 0.12. (8.29) 


This analysis justifies that the soft SUSY breaking induced MSSM inflation is one of 
the examples of inflection point inflation where the present reconstruction technique 
holds good perfectly. 
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9 Multipole scanning of CMB spectra via reconstructed ef¬ 
fective potential 


In this section I study the CMB TT, TE, EE, BB-angular power spectrum 19 . The 
angular power spectra are defined as: 


r ,XY _ 

u e = 


2£ + l 


^ 1 ( a X,im a Y/m) j 


m=—l 


X , Y — T,E,B. 


(9.1) 


Further substituting the inflationary input spectra P{k ) = {Ps(k), Px^k)} and the 
angular power spectra of CMB temperature fluctuations and polarization 


Cf Y = ^ f k 2 dk P{k) A X j{k)A Y e(kl , (9.2) 

Inflation Anisotropies 


where 

/ VO 

dr] S x (k , rj) Pxe(k[r)o - y}\ ■ (9.3) 

Sources Projection 


The integral (9.2) relates the inhomogeneities predicted by inflation, P(k ), to 
the anisotropies observed in the CMB, . The correlations between the different 
A" and Y modes are related by the transfer functions A xe(k) and A Ye(k). The trans¬ 
fer functions may be written as the line-of-sight integral Eq (9.3) which factorizes 
into physical source terms Sx(k,rj ) and geometric projection factors Pxe(k[r ]o - rj\) 
through combinations of Bessel functions. 


19 In this work I have not consider possibility of other cross correlators i.e. TB, EB as there 
is no observational evidence of such contributions in the CMB map. Also till date there is no 
observational evidence for inflationary origin of BB angular power spectrum except from CMB 
lensing. However, for the completeness in this paper I show the theoretical BB angular power 
spectra from the reconstructed potential. 
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Figure 17. We show the variation of CMB TT Angular power spectrum with respect 
to the multipole, l for scalar modes with the choice of best fit reconstructed model pa¬ 
rameters: V(4>o) ~ O(10~ 9 Mp), V'(fo) ~ O(10- 10 Mp), V"(<t> 0 ) ~ O(10- u M 2 ), 
V'"(<fo) ~ 0(1O- 10 M p ) and V""(<f>o) ~ 0(1O -9 ) obtained from Planck 2015 data, which 
is consistent with the bound on held excursion value. From this analysis finally I get: 
P s (k *) ~ 2.215 xlO" 9 , n s (h) ~ 0.962, r(fc*) ~ 0.2, a s (k*) - 10" 2 and n S (K) ~ 5xl0" 3 . 


Further expressing Eq (9.2) in terms of TT, TE, EE, BB correlation I get: 


For scalar : Cj 7 = — j k 2 dk Ps(k) A E e(k)A E e(k) , 


Cg = (An) 2 J k 2 dk P s (k)A Te (k)A Ee (k), 
Cf L = (An) 2 [ k 2 dk P s (k) A E£ (k)A T£ (k ), 


For tensor : Cf = (An) 2 J k 2 dk P E (k) A E g(k)A E e(k) 
Cj T = - f k 2 dk P T (k)A Te (k)A Tl (k), 


n 


Cj E = (4ir f I k 2 dk P T (k) A T ,(t)A H (t), 


Cf E = (4!t)' 2 / k 2 dk P T (k) A E ,(k)A Tt (k ). 


(9.4) 

(9.5) 

(9.6) 

(9.7) 

(9.8) 

(9.9) 
(9.10) 
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Figure 18. We show the variation of CMB TE Angular power spectrum with respect 
to the multipole, l for scalar modes with the choice of best fit reconstructed model pa¬ 
rameters: V((j) o) ~ O(10- 9 Mp), V'fa) ~ O(10- 10 Mp), V"(c /) 0 ) ~ 0{ 10“ n M 2 ), 
V (4>q) ~ O(10 _1 ° Mp) and V (4>o) ~ 0(1O -9 ) obtained from Planck 2015 data, which 
is consistent with the bound on held excursion value. From this analysis finally 1 get: 
P s (k*) ~ 2.215 xlO" 9 , n s {h) ~ 0.962, r(fc*) ~ 0.2, a s {k *)-10" 2 and r, S {K) ~ 5xl0" 3 . 

where the inflationary power spectra {Ps(k), Px(k)} are parametrized at any 
arbitrary momentum scale as: 

/a ns - 1 ^ ln (h)+? ln 2 (s)+-" 

Ps{k) = Ps(K) (-J , (9.11) 

/ b \ n r+^ ln (^) + ^ ln2 (^)+ - 

P T (k) = r(k)P s (k ) = P T (fc*) ( . (9.12) 

It is important to note that the cosmological significance of the E and B decom¬ 
position of CMB polarization carries the following significant features: 

• Scalar (density) perturbations create only d-modes. 
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Figure 19. We show the variation of CMB EE Angular power spectrum with respect 
to the multipole, l for scalar modes with the choice of best fit reconstructed model pa¬ 
rameters: V((j) o) ~ O(10- 9 Mp), V'fa) ~ O(10~ 10 Mp), V"(c/) 0 ) ~ 0{ 10“ n M 2 ), 
V (4>q) ~ 0(1O _1 ° Mp) and V (4>o) ~ 0(1CF 9 ) obtained from Planck 2015 data, which 
is consistent with the bound on held excursion value. From this analysis finally I get: 
P s (k*) ~ 2.215 xlO" 9 , n s {h) ~ 0.962, r(fc*) ~ 0.2, a s {k *)- 10" 2 and r, S {K) ~ 5xl0" 3 . 


• Vector (vorticity) perturbations create mainly 5-modes. However, the contri¬ 
butions of vectors decay with the expansion of the universe and are therefore 
sub-dominant at the epoch of recombination. For this reason I have neglected 
such sub-dominant effects from our analysis. 


• Tensor (gravitational wave) perturbations create both 5-modes and 5-modes. 


To compute the momentum integrals numerically and to analyze the various 
features of CMB angular power spectra from the prescribed reconstruction algorithm 
1 use a numerical code “CAMB” [44], For the numerical analysis 1 use here the best 
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CMB BB (tensor) Angular power spectrum 


CMB TT (tensor) Angular power spectrum 




(a) l(l + 1)C EB / 27 r vs l (tensor) (b) l{l + l)Cf T /2 tt vs l (tensor) 


CM B TE (tensor) Angular power spectrum CM B EE (tensor) Angular power spectrum 




(c) l(l + 1)C EE / 27 t vs l (tensor) 


(d) l(l + 1)C EE / 27 t vs l (tensor) 


Figure 20. We show the variation of (a) CMB BB Angular power spectrum, (b) CMB 
TT Angular power spectrum, (c) CMB TE Angular power spectrum and (d) CMB EE 
Angular power spectrum with respect to the multipole, l for scalar modes with the choice 
of best fit reconstructed model parameters: V(4>o) ~ 0(1O -9 M 3 ), V' (</>o) ~ 0(1O~ 10 iff))), 
V"(<f>o) ~ O(10~ n Mp), V'"{(f> o) ~ O(10 _1 ° M p ) and V""(fo) ~ 0{ 10~ 9 ) obtained from 
Planck 2015 data, which is consistent with the bound on the field excursion value. From this 
analysis finally I get: Ps(£:*) ~ 2.215 X 10 -9 , ns{k+) ~ 0.962, r(/c*) ~ 0.1, asik*) ~ —10 -2 
and ac s(k*) ~ 5 X 10~ 3 . 
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fit reconstructed model parameters: 



v{<h) 

~ £>(10~ 9 Mp), 

(9.13) 

V\<t> o) 

~ C>(10- 10 M 3 ), 

(9.14) 

V”(<M 

~ C>(10- n M 2 ), 

(9.15) 

V"\ct> o) 

~ C>(10~ 10 M p ), 

(9.16) 


~ C>(10“ 9 ) 

(9.17) 


which is compatible with Planck 2015 data. For this specific choice of multi parameter 
space finally the inflationary observables are estimated as: 

Ps(k *) 

~ 2.215 x 10~ 9 , 

(9.18) 

n s {h) 

~ 0.962, 

(9.19) 

r(h) 

~ 0.1, 

(9.20) 

ois{K) 

-10“ 2 , 

(9.21) 

*>s{K) 

~ 5 x 10~ 3 . 

(9.22) 


Additionally I take ACDM background in which I fix: 


f \h 2 = 0.022220, 

(9.23) 

Q c h 2 = 0.119700, 

(9.24) 

f \h 2 = 0.000640, 

(9.25) 

o' 

a 

Cl 

(9.26) 

100 6 = 1.040832, 

(9.27) 

N eff = 3.046. 

(9.28) 

Further I put all these inputs to “CAMB” and modify the inbuilt parameterization of 
power spectrum for scalar and tensor modes accordingly. Finally from the analysis 1 
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get: 


Q m (= 1 -n K - n A ) = 0.314658, 


a8(all matter) = 0.8234. 


rs(zdrag) = 147.35 Mpc, 
^d(z*) = 0.1407 Mpc -1 , 
100 0 D = 0.160737, 

100 6 eq = 0.814740, 

1 Irecomb 280.75 klpC, 

r] now = 14164.5 Mpc, 


= 1089.87, 

rs(z+) = 144.64 Mpc, 
D a (z+) = 13.89228 Mpc, 
Zdrag = 1059.40, 


ZReion 9.959, 

z eq = 3391.45, 


= 0.685342, 


t 0 = 13.814 Gyr, 


(9.29) 

(9.30) 

(9.31) 

(9.32) 

(9.33) 

(9.34) 

(9.35) 

(9.36) 

(9.37) 

(9.38) 

(9.39) 

(9.40) 

(9.41) 

(9.42) 

(9.43) 

(9.44) 


After performing all the numerical computations via “CAMB” finally from our anal¬ 
ysis I have generated all the theoretical CMB angular power spectra from which I 
observed the following significant features:- 

• In Fig. (17), at low i region (2 < l < 49) the contributions from the running 
{ols,olt), and running of running (ks,kt) are very small. Their additional 
contribution to the CMB power spectrum for scalar and tensor modes becomes 
unity (0(1)) within low-/ region and the original power spectrum becomes 
unchanged. As a result the reconstructed model will be well fitted with the 
CMB TT spectrum at low -l region within high cosmic variance as observed by 
Planck except for a few outliers according to the Planck 2013 data release. But 
to update our analysis with the latest Planck 2015 data set I have not shown 
such high cosmic variance explicitly. In case of WMAP9 data for CMB TT 
spectrum the cosmic variance in the low-/ region is not very large compared 
to the Planck low-/ data. But our reconstructed model is also well fitted with 
WMAP9 low-/ data also which I have not shown explicitly in the plot to update 
the analysis using Planck 2015 data. It is also important to mention that if I 
incorporate the uncertainties in the scanning multipole / for the measurement 
of CMB TT spectrum, then also our prescribed analysis is pretty consistent 
with the Planck 2015 data. Further if I move towards high £ regime (47 < 
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I < 2500) the contribution of running and running of running become stronger 
and this will enhance the power spectrum to a permissible value such that it 
will accurately fit high-/ data within very small cosmic variance as observed 
by Planck. In this way one can easily scan over all the multipoles starting 
from low-/ to high-/ using the same momentum dependent parameterization of 
tensor-to-scalar ratio as prescribed in Eq (2.17). 

• From Fig. (17), we see that the Sachs-Wolfe plateau obtained from our pro¬ 
posed reconstructed model is non flat, confirming the appearance of running, 
and running of the running in the spectrum observed for low / region (/ < 47). 
For larger value of the multipole (47 < / < 2500), CMB anisotropy spectrum 
is dominated by the Baryon Acoustic Oscillations (BAO), giving rise to several 
ups and downs in the CMB TT spectrum. In the low / region due to the pres¬ 
ence of very large cosmic variance there may be other pre-inflationary scenarios 
which might be able to describe the TT-power spectrum better. In my study 
I have considered only the possibility for which the behaviour of reconstructed 
model is analyzed for both low and high / regions. 

• From Fig. (17), Fig. (18), Fig. (19), we observe that if I include the uncertainties 
in multipole / as well in the observed CMB angular power spectra then the 
proposed reconstructed model is pretty consistent with the CMB TT, TE, EE 
for scalar mode from Planck 2015 data. 

• In Fig. (20(a)), Fig. (20(b)), Fig. (20(c)) and Fig. (20(d)) I have explicitly 
shown the theoretical CMB BB, TT, TE and EE angular power spectrum from 
tensor mode. Most importantly, if inflationary paradigm is responsible for the 
nearly de-Sitter expansion of the early universe then the CMB BB spectra for 
tensor modes is one of the prime components trough which one can detect the 
contribution for primordial gravitational waves via tensor-to-scalar ratio. But 
till date only the contribution from the lesing B-modes are detected via South 
Pole Telescope [45] and Planck 2015 +BICEP2/Keck Array joint mission [14]. 
But confirming the sole inflationary origin from the detection of the de-lensed 
version of the signal is not sufficient enough to draw any final conclusion 20 . 
There are other possibilities as well through which it is possible to generate 
CMB B-modes, those components are: 

1. Primordial Magnetic Field [47-49], 

2. Gravitational Lensing [45, 50], 

3. CP asymmetry in the lepton sector of particle theory [48, 49] etc. 

20 In this respect one may consider the alternative frameworks of inflationary paradigm as well 

[46]. 
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10 Conclusion 


In this paper I have obtained the following results: 

• We derive the most generalized analytical expressions for the field excursion 
\A(p\/M p in terms of r(k *), 1/(0*), //(&*) for a generic sub-Plackian model of 
inflation by taking into account of higher order slow roll corrections and Buch- 
Davies initial condition. We derive the results for various parameterization of 
the primordial power spectrum by allowing: 

1. scale invariant feature, 

2. further modification in presence of spectral tilt, 

3. the modification in presence of spectral tilt and running of the tilt, 

4. the effect of spectral tilt, running of the tilt and running of the running 
of tilt, 

through which it is possible to evade the well known Lyth bound. Using these 
relations one can easily rule out all classes of super Planckian or trans-Planckian 
inflationary models available in inflationary literature. 

• For the completeness, I also extend the idea by deriving the most generalized 
analytical expressions for the field excursion \ A(j)\/M p in terms of r(fc*), U(0*), 
H(k+) for a generic sub-Plackian model for the various limiting situations of 
non Buch-Davies vacuum. 

• In order to satisfy the observational constraints from Planck 2015 and Planck 
2015+BICEP2/Keck Array joint constraints, one requires a non-monotonic 
evolution of the first slow-roll parameter, ey, if one wants to build a model 
of inflation with a sub-Planckian field excursion and VEV, as pointed out in 
earlier Refs. [21, 23, 24, 26]. 

• Hence in this paper I have reconstructed the inflationary potential around the 
VEV 0 O by computing the Taylor expansion co-efficients U(0 o ),U / ( < / > o);b r// (0o), 
V'"(0o) and V /W (0o) in terms of the Taylor expansion co-efficients V (0*),U / (0*), 
V"(0*), V'"(0*) and U w/ (0*) at the CMB pivot scale 0* using matrix inversion 
technique using Planck 2015 and Planck 2015+BICEP2/Keck Array joint data. 

• In order to satisfy the current observational constraints from Planck 2015 and 

Planck 2015+BICEP2/Keck Array joint data, I have found that the upper 
bound of the scale of inflation to be within: \JV+ < (9(1.856 — 1.926) x 10 16 GeV, 
for the upper bound on field excursion of inflaton varying within: < 

(9(0.223 - 0.242). 
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• We have also derived a new set of second order consistency relationships for 
any generic sub-Planckian model of inflation within slow-roll prescription. In 
particular, this can be treated as a discriminator to break the degenracy be¬ 
tween various cosmological models of inflation. Most importantly, using these 
relations one can break also various class of models of inflation as well. 

• Furthermore, if the data could be refined to constrain tensor spectral tilt tit, 
running and running of the running Kt, then this would play a very crucial 
role to break the degeneracy between various cosmological parameters esmated 
from different class of inflationary models and also rule out various models 
available in inflationary literature. 

• Within our prescribed methodology I have also found that the choice of the 
field interval $ is sensitive to the number of e-foldings AN ~ 0 (8 — 17), which is 
necessarily required to solve the horizon problem associated with standard big 
bang cosmology and to produce sufficient amount of inflation by constraining 
the Taylor expansion co-efficients of the generic potential using Planck 2015 
and Planck 2015+BICEP2/Keck Array joint data. 

• Finally, I have checked the validity of our prescribed reconstruction technique 
by fitting the theoretical CMB angular power spectra from TT, TE, EE for 
scalar mode and BB, TT, TE, EE correlation for tensor mode with the observed 
Planck 2015 data within multipole scanning region 2 < l < 2500, in which I 
cover the momentum scale within the window, 4.488 x 10 -5 Mpc -1 < k < 
0.3 Mpc -1 . 

• Typically, within the framework of particle physics, the nature and shape of 
the inflationary potential will not just be a single monomial. In principle the 
inflationary potential could contain quadratic, cubic and quartic renormalizable 
interactions for effective field theory, or even higher order non-renormalizable 
UV cut-off scale suppressed effective operators arising from integrating out the 
heavy degrees of freedom originated from string theory inspired hidden sector. 
In this respect our derived expressions and numerical results are important for 
reconstructing a particle physics motivated effective potential for inflation in 
a successful fashion. One can extend the prescribed methodology to various 
string theory motivated effective potentials, where the VEV of the inflaton field 
and field excursion both are sub-Planckian. Consequently, the reconstructed 
version of the inflationary potential allow only the concave feature, which is 
compatible with both particle physics and string theory framework. 

The future prospects of the present work are appended below: 

1. First of all, using the non Bunch-Davies initial conditions one can repeat 
the prescribed methology presented in this paper. Then applying the ob- 
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servational constraints from Planck 2015 and Planck 2015+BICEP2/Keck 
Array data one can constrain the exact version of the non Bunch-Davies 
initial condition applicable for effective held theory description of inha¬ 
tionary paradigm. 

2. Hence applying the observational constraints one can derive the hrst order 
and second order inhationary consistency relations from the exact version 
of the non Bunch-Davies initial condition. 

3. Further, one can compare the results obtained from non Bunch-Davies 
and Bunch-Davies vacuum and check which initial condition is more com¬ 
patible with the observed data. This will in turn hx the ambiguity of 
choosing the proper initial condition for inhationary paradigm. Also by 
doing this analysis one can further constrain the nature and shape of the 
inflationary potential. 

4. Using various other parameterization of primordial power spectra one can 
also derive the most generic relationship between tensor-to-scalar ratio, 
scale of inhation, Hubble parameter and the held excursion and also repeat 
the rest of the reconstruction methodology presented in this paper. 

5. As I have already mentioned at the end of last section that one of the 
prime source of generating CMB B-mode polarization is primordial mag¬ 
netic held, by applying the observational constraint one can reconstruct 
an effective held theory of inhationary magnetogenesis for sub-Planckian 
inhation. 

6. Finally using the results for non Bunch-Davies initial condition one can 
also check the compatibility of the theoretical CMB angular power spec¬ 
tra with the temperature anisotropy and polarization data obtained from 
Planck 2015. 
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Appendix 

A. Slow-roll Integration from the reconstructed po¬ 
tential: 

The expressions for the slow-roll parameters (ey, rjy, <7y) can be expressed as: 


ev{k) = ev ~Y ln (f )+ t ln2 (f I + 

Vv(k) = r] V - ^ a ' S o 3ai ' 1 h 1 (+ \ ( K s ~ 3 kt + + «t} 


ns — 3n T — 1 


y (1 - n s ~ 3 n T } ) In 2 ( £ ) + 


(v{k) = & ~ \ («s ~ 4 k t + 4 n 2 T {n s - n T - 1}) In (J^j 

+ \ {16 n% + 9 n T [n s - 3 n T - 1] + [n s - 3 n T - l] 2 + 2^}) In 2 + 


a v(k) = <Tv + <Ty (1 n s ) In ( — 


cr 3 / k \ 

+ -44 (30n| + 20n r [n s - 3 n T - 1] + £y + 2 [n s - 3 n T - l] 2 ) In 2 ( — j + 

( 10 . 1 ) 

where at the pivot scale k = k* the slow roll parameters are defined as follows: 


* - f ( V v ) , w = Ml (U) = Ml (PU) ,4 = Mp 
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Further in terms of the reconstructed potential I get: 


e v 


Vv ~ 

£•2 ^ 
sy ~ 


Ml 

2 V (<j>o ) 2 
M 2 

±v± p 

V(<fio) 2 


V (0o) + V (0o) (0 — 0o) + ^0 _ 0 Q ) 2 _|_ V (4>o) ^ _ 0^3 _|_ 

^ (0o) + (0o) (</> — 0o) + - —(0 — 0o) 2 + • • • , 

^'"(0o)^(0o) + (^""(0o)^'(0o) + ^'"(0 O )^"(0o)) (0 - 0 O ) 


2 


cr 


3 

v 


V(<po) 3 


+ (V"'(<a„)V'(«/.„) + - a,) 2 + •••], 

r (0o) + r (<^ 0 ) (<7 — ^o) + -— — «W 2 + 2 — — ^o) 3 + ■ ■ ■ 


(10.3) 

Furthermore, the inflationary observables, i.e. the amplitude of scalar and tensor 
power spectrum (P$, Pt ), spectral tilt (■ ns,riT ), and tensor-to-scalar ratio (r*) at the 
pivot scale fc* can be expressed as: 


Ps(k*) — [1 _ (2Ce + l)ey + Cg^y] 2 
^(0o) 2 


V 


12tt 2 M61/'(0o) 2 . 


24vr 2 M p 4 ey 
V{<!>o) + ^ (0o)(0* ~ 00 ) + 


,,2V'"M rj . i ,3 


-(</>* — 0o) + 


6 


^"(0o) 

24 


(0* — 00) 4 + ' ' ' 


1 - (2C S + 1) 


v\m 2 mI 

■2V(M 2 


+ Ce 


-( 0 * - 000 


M|V"(0 o) 


v(4>o) 

(10.4) 


Pt(K) = [1 — (Ce + l)ey] 


2 21/ 


3t r 2 M 4 


3tt 2 M 4 


^(0o) + ^ (0o) (0* — 00) + 

, ^""(0o) 


^"(0o)^ ^2,^"(0o)^ J. \3 


-(0* — 000 + 


6 


24 


■(0* — 00) +••• 


3 7/3 , ^v\M 2 Ml 

i - + ij 2y(0o)2 


-( 0 * - 000 


(10.5) 


- 1 ~ (2r/y - 6ey) + • • • 



( 21/" (0o) 

V ^(0o) 


^'(0o) 2 0 0 2V w, (^o) 

270000 0 V ^(0o) 


21/000)l^ ,, (0o) \ 

^(0 O ) 2 ; 


( 0 * 


0o) H- 

( 10 . 6 ) 


nr ~ —2ey + • • • 



70000 


y'(0o) + ^"(0o) (0* 
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00) + • • • 


(10.7) 
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r(h) fa 16e v [1 + 2C E (e v - Vv)} H- 


8Mp 


^ (00 ) + V (00 ) (0* 
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<M + ■ ■ ■ 


1 + 2 C E M 2 p 



V"^ o) Y 

V(0o) ). 
( 10 . 8 ) 


a.s(K) = (jf£^ ~ ( 16 We y - 24e 2 v - 2^) H-, 


8 M„ 4 

^ (0o) + ^ (00 ) (0* — 00 ) + 
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2 r 


V (00) + ^ (00) (0* — 0o) + 


QMp 

WoT 


y'(0o)+y"(0o) (0* 00 ) 


2Mp 

W 
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(10.9) 


«s(fc*) = 
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= i6M 4 ( Y^YlM _ EO ) 4 


+ 
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P V Wo) 3 ^(0 O ) 4 

The crucial integrals of the first and second slow-roll parameters (ey, r/y ) appearing 
in the right hand side of Eq. (2.14), which can be written up to the leading order as: 


d<t>e v & - 


i ” M|+ 2 c p / 0.-0 o y +1 f/ 
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d<f> Vv ~ ^2 
q =o 


Mf 1 B q 

(? + !) 



9+1 
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9+1 


- 1 


(10.13) 

where I have used the (0 — 0 O ) < Mp (including at 0 = 0* and 0 = 0 e ) around 
0o- The leading order dimensionful Planck scale suppressed expansion co-efficients 
(C p ) and (D g ) are given in terms of the model parameters (V (0o), l/(0,o), • • •), which 
determine the liight and shape of the potential in terms of the model parameters as: 


/'(0o) 2 r 21/ ,, (0o)'P'(0o) 2l/(0o ) 3 
/(0o) 2 ’ 1 /(0o) 2 /(0o) 3 ’ 

/"(0o) 2 51/'(0o) 2 ld' / (0 o ) l/ / (0o)/ /// (0o) 

l/(0o) 2 /(0o) 3 l/(0o) 2 ’ 

/ , (0o)/ ,/// (0o) 7f/(0o)V"(0 o ) /"(0o)/ ,/, (0o) 

31/(0o) 2 31/(0o) 3 z(0o) 2 

41/ / (0q)1/' / (0o) 2 91/ , (0q) 3 1/ // (0q) 

/(0o) 3 /(0o) 4 


V\<h) n V"\^) /'(0o)/"(0o) 

/( 00 ) ’ 1 V ( 00 ) Z(0o) 2 ’ 

v""^) V\<t> Q )V"\<f> Q ) /"(0O) Z"(0o)!/(0o) 2 

2/(0o) /(0o) 2 Wo) 2 Wo ) 3 ’ 

41/'(0o)h 2 21/ / (0o) 3 h l/'(0o) 2 l/ / "(0o) 21/ ,/ (0 o )l/ ,// (0o) / // "(0o)l/'(0o) 

Wo) 3 Wo) 4 V (0o) 3 31/(0o) 2 21/(0o) 2 


(10.14) 

Here Wo), /'(0 O ), WW H""(0 o ) ^ 0, 1/"(0 O ) = 0 and Wo),, /'"(0o), V""(0 o ) / 
0, V (0o),H (0o ) = 0 are two limiting situations which signifies the inflection point 
and saddle point inflationary setup. For details, see Ref. [23]. 
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B. Momentum Integration from various parameteri¬ 
zation of tensor-to-scalar ratio: 


In general, the tensor-to-scalar ratio can be parametrized at any arbitrary momentum 
scale as: 


r(k) = < 


r(K 


r (h) ( — 


n r (A:*)—n s (A;*)-(-l 


r (h) ( — 


k xn r (*.)-^W+l+ aT(t,) a 7 (t * 1 H&) 


r (h) ( — 


fc \nr(^)-n 3 (W ar(t * ) a | ag(tt) K £)+ ln 2(_j^ 


for Case I 


for Case II 


for Case III 


for Case IV. 


(10.15) 

where A;* be the pivot scale of momentum. Here these four possibilities are:- 


• Case I stands for a situation where the spectrum is scale invariant. This is the 
similar situation as considered in case of Lyth bound [39]. This possibility also 
surmounts to the Harrison & Zeldovich spectrum, which is completely ruled 
out by Planck-fWMAP9 data within 5cr C.L. 


• Case II stands for a situation where spectrum follows power law feature through 
the spectral tilt ( ns, nr ). This possibility is also tightly constrained by the 
WMAP9 and Planck+WMAP9 data within 2cr C.L. Recently in Ref. [51] the 
authors have explicitly shown that power law feature in the primordial power 
spectrum is ruled out at more than 3cr C.L., 

• Case III signifies a situation where the specrum shows deviation from power 
low in presence of running of the spectral tilt (as, olt) along with logarithmic 
correction in the momentum scale as appearing in the exponent. This possi¬ 
bility is favoured by WMAP9 data and tightly constrained within 2cr window 
by Planck+WMAP9 data, 


• Case IV characterizes a physical situation in which the spectrum is further 
modified compared to the Case III , by allowing running of the running of 
spectral tilt (ks, Kt) along with square of the momentum dependent logarithmic 
correction. This case is satisfied by both WMAP9 and Planck+WMAP9 data 
within 2cr C.L. This is the only criteria which is always satisfied by a general 
class of inflationary potentials. In this article, I have only focused on this 
possibility, from which I have derived all the constraint conditions for a generic 
model of sub-Planckian inflationary potentials. 
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Let us start with the computation of momentum integration where I investigate 
the possibility of four physical situations as mentioned in Eq (10.15) finally leading 
to: 
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for Case I 


for Case II 


for Case III 


for Case IV. 


(10.16) 


where in a realistic physical situation one assumes the pivot scale of momentum 
k * ~ kcmb • Now further substituting Eq (3.7) on Eq (10.16) I get: 
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for Case I 
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(10.17) 


C. Taylor expansion co-efficients of inflationary po¬ 
tential: 


To write down all the Taylor expansion co-efficients in terms of the inflationary 
observables at the pivot scale Ay I start with the slow-roll parameters €y, rjy, £y, <jy 


72 



which can be expressed as: 
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where 


v (</>*) = v(h) = k = -p s (K)mm;, 


x * = 


r/vV 2r (*h) - \ ~as(h)< 


r(K) 


(10.24) 

(10.25) 


Finally we are left with the following Taylor expansion co-effiecients (derivatives of 
the potential) at the pivot scale in terms of the inflationary observables: 
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D. Field excursion from various parameterization of 
tensor-to-scalar ratio: 

Using the result of Eq (10.12), Eq (10.13), Eq (10.17) I get the following expression 
for the held excursion in terms of the tensor-to-scalar ratio and other inhationary 
observables as: 
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Substituting AN ~ (9(8—17) in Eq (10.17) for the previously mentioned four physical 
situations and further using Eq (2.14) the field excursion can be constrained as: 

Planck (2013)+WMAP-9+high L: 


A 0 




< 


(9(0.98 — 2.08) for Case I 
(9(0.87 - 1.88) for Case II 
(9(0.51 - 0.96) for Case III 


(9(0.239 - 0.241) for Case IV . 

Planck (2014)+WMAP-9+high L+BICEP2 (dust): 


A 0 




= 


(9(2.32 - 3.12) 

for Case I 

<9( 1.73 - 2.63) 

for Case II 

(9(0.62 - 0.94) 

for Case III 

(9(0.242 - 0.354) 

for Case IV. 


Planck (2015)+WMAP-9+high L(TT): 
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Planck (2015)+BICEP2/Keck Array: 
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(10.31) 


(10.32) 


(10.33) 


(10.34) 
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In this paper I have only focused on the last possibility, from which I have derived 
all the constraint conditions for a generic model of sub-Planckian inflationary po¬ 
tentials. Also the last possibility is important because within this it is possible to 
generate large value of tensor-to-scalar ratio along with field excursion A fi < M p . 
This also validates the effective field theory prescription within the regime of in¬ 
flationary paradigm. There are other possibilities as well through which one can 
address this crucial issue in the context of inflation. Those possibilities are:- 

• Multi-field inflationary prescription, 

• Randall-Sundrum braneworld [20], 

• Higher curvature gravity [52], 

• Other ghost-free modifications in GR [53], 

• Inflation from torsion [54] etc. 
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